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ABSTRACT 



gjO^T, b . 

The purpose of this investigation was to determine the feasibility 
of using an analog computer to obtain the time derivative of the tempera- 
ture response of a compact heat exchanger surface subjected to a step 
change in incoming fluid temperature; and to investigate the effect of 

the ratio of flow length to hydraulic diameter (L/D ) on the heat 

H 

transfer and flow friction characteristics of compact heat exchanger 
surfaces . 

The method of maximum slope developed by Locke and modified by 

\ t 

Howard to include conduction parameter was used to determine the heat 
transfer information included herein. 

The results show that an analog computer can be a useful tool to 
aid in the collection and reduction of data. Results in the L/D 

H 

investigation were generally inconclusive and bear further investigation. 
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NOMENCLATURE 



English Letter Symbols 

A Matrix total heat transfer area sq ft 

A Matrix minimum free flow area sq ft 

c 

A- Matrix total frontal area sq ft 

fr 

A Matrix solid cross-sectional area sq ft 

Q 

available for thermal conduction 



a 

s 

b 

C 

c 




C 

s 

c 

p 



c 

s 




E 

G 



Fin thickness 

Flow passage perimeter (/3A^) 

Jet contraction-area ratio for 
circular tube 

Fluid stream thermal capacity 
(mc p ) 

Matrix thermal capacity (W c ) 

s s 

Fluid specific heat at constant 
pressure 

Matrix material specific heat 

Flow passage hydraulic diameter 

(4r h ) 

Friction power per unit area 

Flow stream mass velocity (m/A c ) 

Proportionality factor in Newton's 
Second Law 



ft 

ft 

dimensionless 

Btu/(hr deg F) 

Btu/deg F 
Btu/(lbm deg F) 

Btu/(lbm deg F) 
ft 

hp/sq ft 
lbm/(hr sq ft) 

32 . 2(lbm ft) /(lbf sec 2 ) 



h Surface heat transfer coefficient Btu/(hr sq ft deg F) 

for convection; heat transfer 
power per unit area per degree 
temperature difference 

K Contraction loss coefficient for dimensionless 

c 

entrance to heat exchanger 



K 

d 



Momentum velocity-distribution 
coefficient 



dimensionless 



K Expansion loss coefficient for dimensionless 

heat exchanger exit 
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k Fluid thermal conductivity 

k g Matrix thermal conductivity 

L Total matrix flow length 

m Mass flow rate 

P Pressure 

P Matrix porosity (A u ) 

c fr 

q Heat transfer rate 

R Gas constant (53.35 

-air) 

r, Hydraulic radius (A L/A) 

n C 

t Temperature 

u Flow velocity 

V Matrix volume 

m 

W Matrix mass 

s 

x Distance along flow passage 

from the matrix inlet 

Greek Letter Symbols 

Compactness (A/V ) 
tn 

Ratio of orifice diameter to 

pipe diameter (d /d) 
o 

Difference or Change (time, 
temperature, distance, etc.) 

Time 

Fluid viscosity 
Density 

Subscripts 

atm Local atmosphere 

ave Average 

Fluid (gas, air) 



/B 

& 

A 

e 

A 

/° 



Btu/(hr sq ft deg F/ft) 
Btu/(hr sq ft deg F/ft) 
ft 

lbm/hr 
lbf/sq ft 

dimensionless 

Btu/hr 

(ft lbf) / (lbm deg R) 
ft 

deg F 
ft/sec 
cu ft 
lbm 
ft 



sq ft/cu ft 
dimensionless 



sec, hr 
lbm/ hr ft 
lbm/ cu ft 



f 
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i Initial, inlet 

m Matrix, mean 

0 At orifice 

s Solid (Matrix material) , static 

STD Standard (temperature and pressure) 

x Local conditions 

1 Inlet conditions (upstream of matrix and heaters) 

2 Inlet conditions at matrix entrance 

3 Exit conditions at matrix outlet 



Dimensionless Groupings 



Fanning friction factor; ratio of wall shear stress to fluid 
dynamic head. Plotted as a function of Reynolds No. to define 
the surface friction characteristics. 

2/3 °* 

Colburn j -factor (Ng t Np r ). This factor plotted vs. Reynolds 

No. defines the surface heat transfer characteristics. 



N, 



Nu 



Nusselt Number (hD^/k) , a heat transfer modulus 



N. 



Pr 



N, 



N 



St 



N 



tu 



A 

r 



Prandtl Number (/^c^/k), a fluid properties modulus 

Reynolds Number (4r^G ///, ) , a flow modulus 

Stanton Number (h/Gc ), a heat transfer modulus 

P 

Number of heat transfer units (hA/mc ) 

P 

Longitudinal heat conduction parameter for solid 

material (k A /mLc ) 
s s p 

Time parameter (hA9/W c ) 

S S 



11 



ACKNOWLEDGEMENTS 



The author takes this opportunity to express his appreciation to 
Dr. Paul F. Pucci, Professor of Mechanical Engineering, for his patience, 
understanding and guidance in this work. 

The continued support of the sponsor. Naval Ship Systems Command, 
is most gratefully acknowledged. 



12 



1 . 



Introduction. 



One of the methods used to determine the heat transfer characteris- 
tics of compact heat exchangers, that is, those heat exchangers with a 
heat transfer surface area to volume ratio greater than 200, is the single 
blow or transient test technique developed by Locke (12) and further modi- 
fied to include the effects of conduction by Howard (6). This method has 
been used by several previous investigators at the Naval Postgraduate 
School to evaluate the heat transfer and flow friction properties of 
several different materials and geometries. 

In the single-blow transient testing technique, a heat exchanger 
matrix is subjected to a step change in fluid temperature. The response 
of the exit fluid temperature is monitored following this step change. 

The maximum rate of change of the exit fluid temperature is uniquely 

related to the dimensionless number of heat transfer units, N = hA/mc . 

tu p 

Experimentally, this response is recorded as a temperature-time curve and 
the maximum rate of change determined from the maximum slope of this curve 
Location of the maximum slope of this curve has been determined by 
a "cut and try" approach. This method leaves some uncertainty as to the 
exact location of the maximum slope and it was felt that an analog com- 
puter with a differentiating circuit could be used to either determine 
the value of maximum slope directly or indirectly by locating the inflec- 
tion point on the response curve which would be plotted simultaneously 
with the computer output. Therefore, one of the objectives of this in- 
vestigation was to evaluate the use of the analog computer as a differ- 
entiator operating directly on the output of a thermocouple. 

A second objective was to investigate the effect of the flow length 

to dydraulic diameter ratio, L/D , on the flow friction and heat transfer 

u 

characteristics of a compact heat exchanger surface. 
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2. Summary of Theory. 

A. Background 

An analytical solution to the "single blow problem" was first 
presented by Anzelius in 1926. Subsequent work was then done by Nusselt 
in 1927, Hausen in 1927 and 1929 and by Schumann in 1929. (14). Schumann's 

solution to the problem involved an explicit solution of the differential 
equations in the form of two infinite series using Bessel functions. 

These solutions generated a family of theoretical curves which were used 
as a standard. A plot of experimental data was then compared with the 
theoretical solutions, and the theoretical curve best matching the ex- 
perimental data was used to determine the desired information. 

Locke (12), differentiated Schumann's solutions and observed that 
the maximum slope of the exit fluid temperature response curve was a 
unique function of N . 

Howard (6), using a digital computer and finite difference techni- 
ques determined a series of solutions for N vs. maximum slope that in- 
cluded the effects of longitudinal conduction, which had been assumed 
equal to zero in Schumann's solutions. 

B. Theory. 

The following theory treats the problem of determining, for a 
gas flowing through a porous solid, the temperature of the fluid and the 
solid as a function of position and time after the incoming fluid has 
been subjected to a step change in temperature. 

The following assumptions are made: 

1. Fluid properties are independent of temperature. 

2. Fluid flow is steady. 
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3. Homogeneous porous solid. 

4. Thermal conductivities of the solid and the fluid are in- 
finite in the direction normal to flow. 

5. The thermal conductivity of the fluid is zero in the 
direction of flow. 

The following analysis is based on the energy balance shown below: 



™ C P f f [Z=J> h b C - f s )c/x l =£> ihcp (tf + c/x) 



V 

Ciz£> | So<-/p 

MH 



czz[> - k s A,(U* + 




The initial conditions and the boundary conditions for the analysis 
derived from an energy balance on an element of the porous solid are: 

1. The matrix is initially at a uniform temperature. 

2. There is a step change in the incoming fluid temperature 
at time equal to zero. 

3. The matrix boundaries are abiabatic. 

Let: t g = temperature of the solid. 

t ^ = temperature of the fluid. 

The heat rates from the energy balance are: 

1 . Energy absorbed by the solid = A s c s ( 3 ^ ^ d x 

2. Heat transferred to the solid by convection = hb (t f - t s ) cfx 

3. Heat transferred from the fluid by convection = 

™ c p Mr 



4. Heat transferred in the solid by conduction = 
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The energy balances are: 



1. Fluid: (2-1) 

*n\ Cf di~ + h b (tf - b s ) die =. o 
dy 

2. Solid: (2-2) 

fAsC; ^3 dx - | 

Let T' = the generalized time variable (dimensionless) 

T = (hA/W c )(9 - W x/mL) (2-3) 

s s t 

where : 

h = unit conductance for convection heat transfer (Btu/hr sq ft deg F) 

A = matrix total heat transfer area, (sq ft) 

W c = matrix thermal capacity, (Btu/deg F) 
s s 

9 * time, (hrs) 

Wf = mass of fluid in matrix, (Ibm) 
m « mass flow rate of fluid (Ibm/hr) 

x = distance along flow passage in direction of flow measured from 
the inlet, (ft) 

L = total matrix flow length, (ft) 

Multiplying the second term in (2-3) by c /c and rearranging leads to 

P P 

the equation 

T- / lA q _ h_A % / Wf Cp \ 

WsC j -rnLcp l W SCs J 

where the second term may be ignored because W^c is much less than W c 

J b f p s s 

when the working fluid is a gas. Thus the time parameter reduces to 

tA Q (2-4) 

Let z = the generalized position variable (dimensionless) 

z = (hA/mc ) (x/L) = N x/L (2-5) 

P tu 
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and z * N at the matrix exit where x = L. 
tu 

Let X = the longitudinal conduction parameter (dimensionless) 

X = k A /(me L) (2-6) 

s s p 

where: k = thermal conductivity of the matrix (Btu/ hr ft deg F) 
s 

A - Solid matrix cross-sectional area available for thermal 
s 

conduction (sq ft) 

Substituting the above dimensionless groups into the energy balance 
equations (2-1,2) and rearranging, the equations become: 



Fluid: 



Solid: 



Mr - 



hr 



^ N tu i * (*f - t s ) 



(2-7) 

( 2 - 8 ) 



If thermal conduction is assumed zero in the direction parallel to 
flow, equations (2-7,8) simplify to: 



Fluid: 



Solid: 



r 

dz 

iL 






(2-9) 

(2-10) 



Schumann* s solutions to these equations are: 

oo 



fc.~ ti) - /_ e 
(u, - u) 



(z + r) 



7>=/ 



( 2 - 11 ) 



,oo 



j - / - e ’ rz * / ‘ > Z z *s£r 



-n^o 



(2-12) 



Locke (12) then differentiated these solutions for constant N to 

tu 

arrive at the following equation for the slope of the heating curve: 
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(2-13) 





- J, C 2 i 



e 




The response of the downstream fluid temperature is used at x - L 
where z then equals N and t^ equals t^ ^ * 

At low flow rates the effects of longitudinal conduction cannot be 
neglected and these effects were taken into account by Howard (6) who 
applied a finite difference method to equations (2-7) and (2-8) to obtain 
solutions, with the aid of a digital computer, for N as a function of 
both maximum slope and conduction parameter A. . Howard's results were 
plotted and tabulated so that by calculating A and measuring the maxi- 
mum slope of the heating curve, the graph or the tables could be entered 
to determine the corresponding values of N . See Figure 1 and Table 1. 

C. Differentiator. 

Determination of the maximum slope of the exit temperature has 
been done by a visual procedure by previous investigators of this problem. 
Slope was found by sliding a straight edge along the trace of the tempera- 
ture curve until the point where slope appeared to be a maximum was found. 
The tangent line to the curve was drawn at this point and the value of 
slope calculated. It was felt that this method left some doubt as to 
the accuracy of the value of maximum slope and it was decided to use an 
analog computer with a differentiating circuit to determine maximum slope* 
Ideally, the computer output should give the value of slope desired 
directly, but if this method should turn out to be inaccurate, the peak 
of the derivative curve should allow the point of inflection on the 
temperature response curve to be accurately determined, assuming no tim- 
ing errors exist. 
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When using the analog computer to solve differential equations, the 
computer is normally programmed to solve equations by integration, avoid- 
ing differentiation whenever possible because of the inherent noise pro- 
blem. If noise is represented by a sinusoidal such as Asin (£Jt), 
where equals the angular velocity, then its derivative becomes A OJ 
cos ( CJ t) and the magnitude of the differentiated noise signal is direct- 
ly proportional to frequency. With the output of the thermocouples used to 
monitor temperature response in the vicinity of five millivolts and less, 
high frequency noise signals could easily obliterate the data signal 
after differentiation. 

Since the signal in the circuit was D.C., and no alternating cur- 
rent was present, it was felt that the noise problem was minimal and an 
analog differentiator could be used. 

Several circuits are available for analog differentiation (5), the 
simplest being the inverse of the integrating circuit: 

,,c 



e-. o- 





operational amplifier 

Integrating Circuit 



R 




OPERATIONAL amplifier 
Differentiating Circuit 



By summing currents into the amplifier (current through the amplifier 
assumed zero) the following differential equation describes the circuit: 
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c 

^ dt 



or 




dle^ p 

dL"t Rc 



O 



where if we let e^ be the output of the downstream thermocouples, then 

the output of the circuit (e^ above) is equal to the derivative, modified 

by the coefficient (-RC). Taking the Laplace transform and putting the 

equation into the s domain, where 

s = a complex variable 

E (s) = - RC s E . (s) . 
o i 

Letting the product of RC equal to one reduces the equation to: 

E (s) = - s E. (s) . 
o i 

Recalling that the transfer function E o (s)/E i (s) = s is the deriva- 
tive, the equation for the circuit can be rearranged to give: 



Eo (J) 
E ; (s') 



-RC s 



or 



E b (a _ 

El (s) 



- S 



when the RC product equals unity. This circuit is the true derivative 
and is not satisfactory for use if any noise is present, due to noise 
amplification. 

Additional circuits that attempt to overcome the noise problem are 
combinations of differentiators and first order, low pass filters. An 
example of such a circuit is the basic differentiator with additional 
resistors or capacitors, shown below with their corresponding transfer 
function: 



ft, 

e,o — Mr 







R 




« e 0 



_ _ RzCs 
E; Cs) & t Cs+l 



- _ RC. .c 

Ei ^>' Rc 2 s + 1 
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A drawback of the one amplifier circuits shown above is the necessity 
of knowing the frequency of the noise and precalculating the required 
resistor and capacitor values based on eliminating this noise. 

The modifications of the basic differentiator are referred to as 
approximate differentiators. Several of the circuits were tried in the 
course of this investigation and were found unsuitable. One circuit that 
was satisfactory and is recommended for general use (5), can be found in 
Jackson (7). This circuit combines an adjustable low pass filter with 
the differentiator and is as follows: 




The differential equation that follows describes the problem to be 
solved, that is, an approximation to the derivative. It is: 

(t- a) 2 ^ * 2 = (2-14) 

al& Q ft 

If (a) is allowed to approach one, the term containing dz/dt diminishes 
and the equation becomes z = dx/dt when (a) is equal to one. Integrat- 
ing equation (2-14) results in an equation that describes the analog cir- 
cuit shown above: 
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(/- a) z / % dt = x 



(2-15) 



or, rearranging 



Z= X^Q2 -\Zdt 



(2-15a) 



Taking the Laplace transformation and putting (2- 15a) into the s 
domain results in: 



Z (s > * X ^ + CL 2 (S) - Z <0 



(2-16) 






(2-l6a) 



The transfer function for this circuit is: 



Z to _ ( 4 ) = ^ 

X(i) £"/ ft) 

which reduces to: 



E°_(j) __ 
(s) ~~ 



(2-17) 



(2-18) 



when (a) is equal to one. The features of this circuit which add to its 
convenience are the potentiometer, which allows adjustment of the low 
pass filter to filter out the desired frequencies without precalculation 
and the sign changer, which provides the derivative with the proper sign. 



D. Effect of L/D . 

H. 

There is little information published on the effect of L/D„. 

H 

Analytical solutions exist for simple geometries such as circular tubes, 
annuli and parallel planes for the laminar flow fully developed tempera- 
ture and velocity profiles. For complex geometries however, the 
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differential equations become too difficult to solve. From circular 

tube theory it is known that the hydrodynamic and thermal entry lengths 

are a function of Reynolds No. and x/D, where x is the dimension in the 

direction of flow measured from the tube entrance. Kays and London (9) 

list solutions for fully developed velocity and temperature profiles for 

simple geometries that hold for L/D greater than 100. For fully devel- 

n 

oped flow, the ratio of the Stanton No. , Prandtl No. product to the fric- 
tion factor is a constant. This means that on a plot of f and j versus 
N d the two factors should plot as a pair of parallel lines. Since fric- 
tion factor is inversely proportional to Reynolds No. in the laminar 
region, for an arbitrary long tube, friction factor approaches a slope 
of -1.0., and the corresponding j - factor does also. 

Kays (9) has a plot for triangular flow passages where he shows an 
analytical solution for an isosceles triangular shaped passage with L/D 

H 

at infinity. Here the slopes are equal to minus one. 

By using several cores of the same material and flow passage dimen- 
sions and geometries, with only the length and subsequently L/D vary- 
ing, one should be able to experimentally determine the ratio of L/D 

H 

where slopes of the f and j curves are equal and where the flow becomes 
developed to the point that the exit and entrance effects are negligible 
in comparison with the effects of the internal passage. 
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3. Experimental Techniques. 

The experimental apparatus was designed to conform to the idealiza- 
tions stated by Howard (6) so that Howard's conduction parameter curves 
could be used. Briefly, these idealizations are: 

1. Fluid flow in the matrix is steady with a uniform velocity 

and temperature profile at any cross-section. Matrix thermal conductivity 
is infinite normal to the flow direction and finite in the direction of 
flow, making the problem one dimensional in space. These idealizations 
have been met through a specially designed entrance nozzle, flow straight- 
ening screens and the heater wires distributed evenly across the flow 
channel. Conductivity is accounted for in ^ , the conduction para- 

meter. 

2. Large matrix thermal capacity in comparison with that of the 
contained fluid. Use of a gas as the working fluid satisfies this ideal- 
ization and eliminates time dependent terms from the equation for the 
fluid. 

3. Constant and uniform thermal properties of the fluid and 
the matrix. 

4. The convective heat transfer coefficient is some suitable 
average and remains uniform and constant. The idealizations of steps 
three and four are met by restricting the temperature changes to small 
values (about 20 deg F above ambient) where the variation in properties 
can be neglected. 

5. A step change in fluid temperature is imposed at time 
equal to zero, after the matrix and entrained fluid have reached a 
steady state temperature. The heaters are of 0.0031 inch diameter nich- 
rome wire with a response time of less than ten percent of the matrix 
response, giving essentially a step change. Heater time constant for 
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a flow rate of 1,000 lbra/hr is 0,0425 seconds compared to a 1.4 second 
time constant for the matrix response. At low flow rates, about 25 lbm/ 
hr, the heater time constant is about 0.3 seconds which compares to a 
time constant of 57 seconds for the system response. Therefore, the 
heaters provide a very close approximation to a step function. 

A schematic diagram of the experimental system is shown in Figure 3, 
which indicates the location of instrumentation. Figures 4 and 5 are 
photographs of the system. Flow is induced into the system through an 
inlet bell, followed by a flow straightening section, the heaters, the 
matrix test section, and a flow measurement device. Pressure taps mea- 
sure the static pressure at the inlet to the test section, pressure dif- 
ferential across the test section, static pressure in the pipe before the 
flow metering orifice and pressure differential across the orifice. 
Thermocouples are located in the inlet to the apparatus, between the 
heaters and the matrix, immediately downstream of the matrix, and in 
the pipe preceding the orifice. 

Pressure is measured by draft gage for low flow rates and water mano- 
meters for the flow rates over about 350 Ibm/hr. Atmospheric pressure is 
measured on a standard mercury barometer. 

The temperature response recorded is the difference in temperature 
between the inlet temperature, t^, and the matrix exit temperature, t^. 

A more complete description of equipment is found in Appendix A. 

A test run is accomplished by predetermining desired flow rate, cor- 
responding approximate pressure drops and required heaters. Air is drawn 
through the apparatus, flow is adjusted until the desired pressure drop 
across the orifice is attained and the heaters are energized. The matrix 
and the heated air are allowed to come to steady state temperature, then 
power to the heaters is secured and a recording of (t^-t^) as a function 
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of time is made on one channel of a dual channel strip chart recorder. 

The thermocouple output is fed to a differential amplifier, amplified by 
1000:1, then fed to the differentiator on the analog computer. The dif- 
ferentiated signal is then led to the second channel on the recorder where 
it is simultaneously plotted with the time temperature response curve. 

The following information is recorded for each run: inlet static 

pressure, matrix pressure drop, orifice static pressure, orifice pres- 
sure drop, atmospheric pressure, air temperature at the orifice, diameter 
of the orifice plate, ratio of orifice diameter to pipe diameter, chart 
speed of the recorder and corresponding scale factors for the recorder. 
After completion of all desired runs, the values to compute slope are 
taken from the recording traces and included on the data sheet. The 
data sheet layout conforms to the data input section of a digital comput- 
er program that reduces all data and calculates the desired information. 

The derivative was calibrated in the following manner: First, a 

known potential from a potentiometer was used to calibrate the channel of 
the strip chart recorder used to record the derivative. This calibra- 
tion was made in the range of 100 millivolts to 5 volts; which was the 
range of the output of the differentiator on the analog. Again, using 
a known potential, the amplifier was adjusted to give exactly 1000:1 
amplification. A ramp function was set up on the analog, and the output 
of the ramp recorded by the calibrated channel on the recorder. An ac- 
curate measurement was made of the slope of the ramp function from the 
chart trace. The ramp function output was then fed to the differentiat- 
ing circuit and the result recorded on the recorder. From equation (2-14) 
or (2-17) it can be seen that if the potentiometer in the differentiat- 
ing circuit is set equal to one, the circuit reduces exactly to the 
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derivative. However, this was not the case in actuality. 

It was found that as (a) approached one the circuit became more 
unstable, showing up as large oscillations and high overshoot about the 
value of the derivative, and when equal to one, all stability was lost 
and no value for the derivative could be determined. By using the pre- 
viously measured value of slope for a standard, the potentiometer was 
adjusted as close to one as possible, about 0.95, while holding oscil- 
lation and overshoot to a reasonable amount. This resulted in a final 
value of the derivative equal to 99.27, of the value used as a standard. 

Noise is a problem in any derivative circuit and the circuitry used 
for this investigation was no exception. The heating system used in the 
past provided current to the heaters from an A.C. source with voltage control- 
led by a Variac autotransformer. The thermocouples used in the test ap- 
paratus are non-shielded and no shielding was provided on the leads from 
a common terminal strip to the outlets located on the base of the equip- 
ment. In addition, no system of grounding was used in the thermocouple 
circuit. The Moseley dual channel, strip chart recorder used has the 
circuitry to filter out spurious noise signals when this noise was imposed 
on the output thermocouple signal, but when this signal was differentiat- 
ed, the noise level was great enough to completely mask the derivative. 

It was then discovered that much of the interference stemmed from in- 
duced 60 cycle noise from the heater circuitry which is in close proxi- 
mity to the thermocouple circuitry. This was eliminated by using D.C. 
for the source of power to the heaters. There is still considerable noise 
in the system, but the filter on the Moseley recorder is satisfactory to 
handle most of it. For example, on a run where the derivative measured 
178 millivolts/second, the measured noise was 3 millivolts/second peak 
to peak. 
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Amplification of the thermocouple response was necessary so that 
the derivative trace could be plotted, especially at low flow rates. 

At low flow rates the derivative is often less than 0.1 mv/sec which is 
only 10 percent of the lowest scale on the recorder. In order to boost 
the analog output to the more usable ranges of the recorder, an amplifica- 
tion of 1000:1 was used. 

With the amplifier and the analog computer circuits paralleling the 
direct input to the recorder, and because both circuits used the same 
source voltage from the thermocouples, which was only three millivolts 
it was felt that possibly there may have been some interaction between 
systems. Several comparison runs were made to see if this situation did, 
in fact exist. Data was collected for the same flow rate for the follow- 
ing conditions: 

1. Amplifier and analog computer in the circuit. 

2. All components separated from the circuit and direct 
current used for heating alone. 

3. All components separated from the circuit with alternating 
current used for power to the heaters, which duplicates 
the method of previous experimenters. 

All methods of determining maximum slope yielded values of 
within the range of experimental accuracy. 
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4. Description of Test Matrices. 

The heat transfer surfaces used for this project were all of a 
similar geometry and construction, consisting of triangular or corrugat- 
ed fins made from solid sheet stock, and splitter plates made from the 
same material. 

The matrices used during the evaluation of the analog technique had 
been evaluated previously by Ball (2) and Bannon (3). These were used to 
provide a correlation on results and to insure proper experimental techni- 
que was used, that is, that reproducible data was obtained. 

For the investigation of L/D effects, a matrix with a small hydrau- 

ri 

lie diameter was required so that a large value of L/D could be obtain- 
ed within the limiting dimensions of the test apparatus. Fortunately 
a crimping roller that produced a forty fin per inch by eighteen mil 
high triangular fin was available. One mil brass shim stock was used to 
construct the fins and splitter plates. Test cores were made to the 
following lengths: 1/2 in, 3/4 in. 1.0 in., 1% in., 2 \ in., and 3.0 in. 

A stainless steel reference matrix which had been constructed with the 
same equipment and had previously been tested was available for use as 
a comparison standard for flow friction results. 

Further information on core geometries and properties is shown in 
Figures 6 through 10. 
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5. Presentation of Results. 

For each matrix tested, the heat transfer and flow friction char- 
acteristics have been computed, plotted and tabulated. Computed results 
are shown tabulated in Tables II through XI, and are plotted as Colburn- 
j vs Reynolds Number, and Fanning friction factor vs. Reynolds Number. 
Reynolds Number was calculated on the basis of hydraulic diameter for 
each matrix. 

For the investigation of L/D effect, six cores were tested. Fig- 

rl 

ures 15 through 20 show individual data from each core and Figures 21 and 

22 are compilations of f and j data respectively from all the cores. 

Three cores tested by previous investigators in this project. Solar 

No. 1, Solar No. 4 and the stainless steel plate-fin reference were used 

for a basis of comparison of the analog technique, A plot of "f" and 

"j" vs. N for these cores is shown in Figures 11, 12 and 14 where the 
R 

results obtained by the use of the analog are compared with previously 
established information. 

During the investigation it was felt that the response of the 
Moseley recorder was not sufficiently rapid at higher flow rates to pro- 
duce accurate results. A comparison was made between the Moseley record- 
er and a Brush recorder using a simultaneous recording technique for the 
testing of one core sample. The results of this comparison are shown in 
Figure 24. 

One new core was evaluated during the course of this investigation. 
Designated Solar No. 6, this core was a stainless steel, triangular fin, 
plate-fin construction with a small hydraulic diameter (0.00126 ft) and 
a fin height of 0.022 inches. Data for this core is shown in Figure 13. 
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6 . 



Discussion of Results. 



For the determination of the practicability of using the analog 
computer, several cores previously investigated were chosen to act as 
standards of comparison. These were Solar No. 1, No. 4 and a stainless 
steel reference matrix. Figures 11, 12 and 14 show f and j data versus 
Reynolds Number for both methods of determining maximum slope: direct 

reading of analog computer output and the plotting of the tangent line 
at the point of inflection. Plotted on these figures are the data of 
previous investigators at this facility. It was hoped that this data 
could be duplicated so that an accurate comparison could be made. In 
addition, in the L/D investigation, the data presented shows the results 
of both methods used. 

Examination of Figures 11, 12 and 14 show that flow friction results 

matched exactly and that the heat transfer data match in the region of 

higher N , where N is generally above 75. There is a disagreement in 
K R 

the lower flow rate region where the data presented ceases to follow the 

expected straight line behavior, but reaches a maximum and then decreases. 

However, closer examination shows that the information presented by both 

methods of determining maximum slope do agree and in fact, several points 

actually coincide. Inspection of the data traces while reducing data 

showed that the cores with a larger hydraulic radius than that used in 

the cores for the brass L/D investigation, namely Solar No* 1 and No. 4 

H 

indicate two possible maximum slope values when the flow rates are over 
about 500 lbm/hr. The trace of the derivative of the response curve shows 
an initial peaking at time equal to 0+. This curve then goes through a 
minimum point and goes through another peak before it decays to zero. 

If slope is computed at both of these peaks and the corresponding 
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j-factor plotted, it is seen that the values at the later peak follow 
the trend of the earlier data, where no ambiguity exists, while the 
initial peak gives erroneous results. Kohlmayr (11) has extended the 
maximum slope method to include the effects of deviation from the ideal- 
ization of a step input. In his paper he describes the double peaking 
of the derivative of the generalized heating curve, and contends that 
the later peak is the correct value of maximum slope to use. The simul- 
taneous plotting of the derivative gives an immediate indication of the 
location of the exact maximum slope and in this case is an aid in deter- 
mining the proper value of the heat transfer coefficient. 

When flow rates increase, the two peaks of the derivative curve move 
closer together until they are indistinguishable. It was felt that the 
recorder in use had a response time too slow to distinguish between these 
values; in order to check out this contention, a simultaneous recording 
was made on the Brush recorder, which has a more rapid response time 
than the Moseley recorder, 0.0035 seconds to 0.23 seconds respectively 
for rise time. No discernable difference between traces could be noted. 
Actually, at the start of the investigation it was thought that the 
initial peaking was due to a transient response that exceeded the capa- 
city of the recorder to handle, it was felt that the use of the Brush 
recorder would differentiate between signals. It was later in the course 
of this investigation that Kohlmayr'.s results became known. It is in- 
teresting to note that at the upper flow rates, about 1000 lbm per hr., 
the thermocouple time constant is 0.43 seconds compared with a time 
constant for the downstream response of 1.4 seconds. This could possi- 
bly be a source of error at high flow rates. 

The investigation of the effect of L/D on flow friction and heat 
transfer did not produce as good a result as hoped for. Above a value 
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of L/D , hopefully to be determined, fully developed flow exltt* amid 
the plot of j and f data should be two parallel lines with the same slopes. 
In other words, j/f should be a constant. For friction, since f is in- 
versely proportional to N , the slope of the f curve should approach minus 

one as the limiting value as L/D„ becomes arbitrarily large. The slopes 

H. 

of the f and j curves were calculated and compared for each value of L/D . 

n 

It can be seen that instead of approaching a limiting value the slopes 
went through a maximum with the 1„5 inch core, then decreased somewhat. 

All cores were made by this experimenter except the 3.0 inch and the 1.5 
inch cores. Each core was constructed of the same number of fin and 
splitter plates, and the dimensions of the individual pieces were con- 
stant except for flow length, which varied according to the size core 
under construction. In spite of these precautions, the overall stacked 
height of each core varied sufficiently to affect the geometric constants. 
With only length changing and all other factors identical, porosity, com- 
pactness, frontal area, free flow area and conduction area should have 
been identical for all cores, but all varied. 

The core material, brass shim stock, was prepared by slicing on a 
paper cutter. This produced burring of the cut edges, which in turn af- 
fected the flow characteristics. The previously constructed cores did 
not have burred edges. 

Brass was used primarily because a quantity was on hand and the addi- 
tional amount needed was easily obtainable from local suppliers. It was 
possible to construct and test cores from material on hand while waiting 
for additional materials to arrive. Unfortunately, brass has a high 
thermal conductivity which gives a' high value to conduction parameter at 
low flow rates. In addition, conduction parameter is inversely proportion- 
al to flow length and flow rate. These three constituents combine to give 
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high values of conduction parameter at low flow rates and with low flow 

length cores. An inspection of Howard's solutions (Figure 1, Table I) 

show few solutions in the high ,A region, that is, .> 0.1. 

With high flow rates, maximum slope decreased and the solutions gave 

values of N in the range of 1.0 to 3.0, where the possible error is 
tu 

greatest (Figure 2). 

While these problems do not account for the "hump" noticed in the 
j-curve, it does give possible reasons for the spread of data, especial- 
ly for the shorter length cores. In retrospect, it appears that a study 

of L/D„ effects could better be done by using a low r, core to obtain 
H h 

data for larger L/D ratios and a dimensionally similar but with a 

H 

larger r, core to obtain data for small ratios of L/D ; preferably 
n H 

constructed of a low conductivity material. 
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7. Experimental Uncertainties. 



There are several sources of possible error which arise from not 
exactly meeting the idealizations and boundary conditions set forth 
previously. These possible error sources are difficult to accurately 
assess a numerical value to, and will not be discussed further in this 
section. 

Error sources which can be evaluated include uncertainties in 
physical constants, inaccuracies in the determination of geometrical 
constants and inaccuracies in instrumentation. Single runs were usually 
made for each data point except where the validity of a run was in ques- 
tion and a check run was made. This system of taking data does not allow 
a large store of statistical data to be collected where inference to error 
may be made from statistical techniques. In view of these limitations, 
the technique for determination of uncertainties developed by Kline and 
McClintock (10) was used. 

The physical constants for the core materials tested were previously 
determined by Ball (2) and Bannon (3) and the uncertainties listed in 
those references are: 

Physical constants 

k + .5% 
s — 

c + . 57. 
s — 

c + .57. 

P ~ 

N pr + 2.07. 
yU + 1.07. 

Similarly, uncertainties in dimensions have been determined to be: 

A, k c , A , A +1.07. 
fr c s — 

L + 0.57. 

W < 0.17. 

s 
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The weight of the matrix was measured on an analytical balance and any 
uncertainty can be considered negligible compared to other physical 
dimensions and constants. 

Inaccuracies in instrumentation were generally established to be of 
the order of one half of the smallest division of the scale on the in- 
struemnt in use. The largest source of inaccuracies were in the mano- 
meters used for pressure measurement. All temperature recordings for 
conditions at the matrix were differences in temperature and were record- 
ed in inches on the strip chart recorder. The measurement of the tempera 
ture at the orifice was obtained from the reading of the thermocouple 
output on a Rubicon potentiometer. The estimated error of the potentio- 
meter was estimated as approximately 0.1 deg F which reduces to about 
0.2% uncertainty. 

The range of pressure measurements was such that several different 
measuring instruments were required, from draft gages at low flow rates 
to a 120 inch water manometer at high flow rates. Because of this range, 
the maximum uncertainty among all readings was used for the analysis. 
These uncertainties are listed below: 

P + 1 . 25% 
o — 

A p „ ± 

A P +1.70% 

m — 

P _ + 0.005"Hg (negligible) 
atm — 

Uncertainty in measuring maximum slope is estimated to be 2.0% (14). 
This uncertainty in maximum slope is then considered with N and ^ 
to determine the error in N from Figure 2. Close examination of Fig- 
ure 2 shows a large uncertainty in N for large values of conduction 
parameter at high N and for all values of conduction parameter when 
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N is near 2.0. 
tu 

Two examples of uncertainty in Colburn j -factor have been cal- 
culated, one for high N and one for low N D (14). Using Figure 2 and 
the methods of reference (10), the following uncertainties in Colburn j 
have been determined: 

N =3.0, A = 0 N error = + 7.2% j error = + 7.5% 

N tu = 25.0, = 0.05 N error * + 10.0% j error = + 10.2% 

By similar analysis, uncertainties for m, N and f were determined 

K 

to be 1.0%, 2.3% and 4.3% respectively. 
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8 . 



Conclusions . 



The use of a differentiating circuit on an analog computer can be 
used to obtain the maximum slope of the generalized heating curve of the 
transient testing technique for evaluating heat transfer properties of 
compact heat exchanger surfaces. The computer output can be used direct- 
ly to determine maximum slope and indirectly to locate the correct point 
of inflection when multiple maximum slopes occur. 

The effects of L/D^ bear further investigation. The results obtain- 
ed in this investigation were inconclusive and at best represent heat 
transfer and friction characteristics of the individual cores tested. 
Better control of geometry and construction is needed to eliminate error 
caused by geometrical inconsistencies. 
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9. Recommendations for Further Study. 

An investigation should be made in the low N region to determine 
whether the maximum in the j-data is actually there or if an error in 
technique exists. 

Further investigation of L/D„ effects can be made. It would be 

n 

recommended to obtain suitable core samples from an established manu- 
facturer to insure conformity between samples. For investigations in 
the low L/D range, a core with a similar geometry but larger r should 

n u 

be obtained as the small r^ brass cores used were too easily affected 
by high flow rates. 

The cyclic test technique developed by Traister (17), should be 

used for low L/D u cores where N is in the maximum error range, 
n tU 

Further extension of Howard's conduction parameter curves in the 
high N tu » high lambda region is needed. Insufficient data exists in the 
region where lambda is between 0.1 and 10.0 and is at high N . The 
subroutine which interpolates in this region in the data reduction com- 
puter program often cannot handle these values and gives incorrect re- 
sults . 

A system of shielding and grounding should be accomplished along 
with installation of thermocouples which have a response time less than 
that of the Nichrome heaters. This would make the heaters the control- 
ling factor in the response to a step function rather than the thermo- 
couples which at present cannot handle higher flow rates accurately. 

It is suspected that the attenuation noted by Traister in the cyclic 
technique at high flow rates and increased frequency may be due in part 
to the slow response of the thermocouples. Shielding the thermocouples 
and installing a ground system would help alleviate the electronic noise. 
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especially that produced by the linear accelerator located directly 
underneath the test apparatus. 
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<ScM EMATIC Of Ti4T Apparatus 
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Idgure 5. Pest Jection Jbowin Inlet ^ell, Heaters -.nl Latrix Holder. 
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SOLAR 1.0.1 




Latrix I^aterial 




Nickel 


Specific Heat, 


(Btu/lbm deg P) 


0.106 


Thermal Conductivity, k g 


(dtu/hr ft deg P) 


36.0 


I .ate rial ihickness 


(inches) 


0.005 


Heat r-n3fer Area, v 


Uq ft) 


12.768 


Flow Length, L 


(feet) 


0.24917 


Frontal vrea , 


(sq ftj 


0.07022 


Conduction Area, A 

3 


Uq ft) 


0.01078 


Free Flow Area, A 

c 


Uq ft) 


0.05944 


/olune, \[ 


(cu ft; 


o.OI 75 


Compactness, 


(sq ft/cu ft) 


729.7026 


Porosity, p 




0.8464 


Hydraulic Ji one ter, J 


(feet) 


0.0046395 


/eight, / g 


(lom) 


1.5072 



Figure 6. Physical and lecmeirie irroperties cf .iolar 1. 
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-)Ol: lit uO • 4 




- . xtrix Late rial 




430 Jtainle 


-Specific :: -at , c 


(-' w 0.-13 .} 0 j\ 


a. 11 


dermal Jonductivity , k ^ 


( ft ie f) 


1 2.8 


material Tli ickneao 


(inc .es ) 


0 • o kj ^ 


j'l077 L ; ;th, L 


(fc t; 


a. 2441 7 


frontal Area, .V. 

I r 


( oq ft) 


P.C2713 


Volume, / 


(cu ft) 


c.UOS 


Joniuction Vr *a, A 


(oc ft) 


p .01 229 


i*Yee / low ire a, A 

c 


( ooq ft) 


e .05464 


Keat iranafer i.r;a, A 


(3Q. ft) 


14.2765 


I’orojitv, > 




v.L 1 ,<;t 


JO' % pacfcne 


(aq ft/au ft) 


7C.9;2 


^ydraulie Jiameter, J 


(f at) 


v • /5 


,/ei lit, / 


(l:n) 


1 .45216 



i?i rare 7 . i . . 



ical ail-* 



etric i roj'Oi 



O' 



4 . 
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. .atrix } aterial 
Jpecific .ieat, c 

a 

Thermal Conductivity, 

- aterial rhickness 

Flow i.enjth, L 

Front xl Area, A „ 
fr 

/olune, ! 

Free Flo 7/ Area, A 

c 

Conduction Area, A o 
’.eat Jransfer \rea, a 
loroaity, p 
Compactness, 

Hydraulic diameter, J 
•lei lit, .1 

3 

Fiyure 8. Icysical an- 



(jtu/ltra de^ F) 

k ( tu/ hr ft de^ F) 
s 

(incho ;) 

(feet) 

(sq ft) 

(cu ft) 

(sq ft) 

(jq ft) 

(sq ft; 

(sq ft/cu ft; 

(fe t) 

( 1 Dm ; 

Jeoi etric i roperties of 



430 Jtainless 
0.11 
12.8 
0.002 
0.2398 
0.065208 
0.015622 
0.05081 
0.01439 
38.442 
0.779 
2465.0 
0.oC1262 
1 .421 98 

jolar i o. 6. 
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iiV 


r- - , . . . » » ' ' ■ 1 1 

i 1 . ij ; j .* x*. Vj Jv*sy.i 


a taX 


, atrix 1 ate rial 




320 itainless 


Jpecific eat, 


( tu/i r d 


0.11 


Vernal Jor:iuctivity , 


( tu/ - f. ae S) 


7.0 


i. aterial JhicAre >3 


(inc .eo / 


0 . 1 


Slow leu ;th, x. 


U~et, 


o . 1 6667 


Itrontal Area, 


(jq ft. 


C. 070 23 


Volune, 7 


(cu ft; 


o. 01170 


Jonduction ^rea, 


. tcj 


o.^/p2 


tree Pier/ Area, a 


(sq ft) 


o.o94S4 


c 

Heat transfer \rea, A 


(oq it) 


3C . 1 20 


I oro >ity, p 




0.393 


Jonpactnea 


( j i ft/cu ft) 


- 575.t 


Hydraulic diameter, J 


(feet) 


^ . oO 1 3t<3 


,/eifjht, 7 


d'-') 


0.51 93 



ji'jure 9. rhysical and Jeo '.etric Properties of S o jfcairu. 

teference itrix 
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BJ1A3S L => 1.5" 




Matrix Material 




70-30 Brass 


.Specific Heat, c 

s 


(3tu/lbm deg F) 


0.092 


Thermal Conductivity, 


(Btu/hr ft deg F) 


0.092 


Material Thickness 


( inches ) 


C.001 


Flow Length, L 


(feet) 


0.125 


Frontal Area, 

fr 


(sq ft) 


0.07111 


Volume* V 


(cu ft) 


0.CC8QS8 


Free Flow Area, A 

c 


(sq ft) 


0.063534 


Conduction Area, A 

s 


(sq ft) 


0.007577 


Heat Transfer Area, A 


(sq ft) 


22.221 


Porosity, p 




0.89345 


Compactness, 


(sq ft/cu ft) 


2500.15 


Hydraulic Diameter, D 

II 


(feet) 


0.0013998 


/eight, 7 

s 


(lba) 


0.4956 


Fi;:ure 10* Physical and Geometric Properties of Brass 


Gore 1 = 1.5" 
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SOLAR liO. 1 
nColbum-j Factor and Panning Friction 
"Factor vs Reynolds Number 
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Figure 11. Heat Transfer and Flow Friction Characteristics of Solar 1. 
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jj* 30LAR NO. 4 

Colburn- j Factor and Fanning Friction 
Factor vs Reynolds Number 




0.00 
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Figure 1 2. Heat Transfer and Flow Friction Characteristics of Solar 4. 
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Figure 15. Heat Transfer and Plow Friction Characteristics of Solar 6. 
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'ijure 14. Heat Transfer and Plow Friction Characteristics of 
Stainless Steel Plate-Pin Reference katrix. 
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figure 15. Heat Transfer and Plow Friction Characteristics of 
Brass Core, L= 0.5" 
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BRASS L= 0.75" 

Colburn- j Factor and Fanning Friction 
Factor vs Reynolda Number 
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Figure 16. Heat Transfer and Flow Friction Characteristics of 
Brass Core, L =0.75 "• 
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BHA33 L=t.O* 
Oolbuta-j Factor and 
Factor, tb Reynolds Humber 
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Figure 17* Heat Transfer and Plow Friction Characteristics of 
Brass Core, L = 1 *0 " 
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Figure 18. Heat Transfer and Flow Friction Characteristics of 
:Sras3 Core, L =1.5" 



61 





1.0 



0.1 



0.01 



0.001 




10 



100 



1000 



N 



R 



Figure 19. Heat Transfer and Plow Friction Characteristics of 
Brass Core, L =2.25" 
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Figure 20. Heat Transfer and Flow Friction Characteristics of 
Brass Core, L= 3.0" 
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Figure 21. Combined Heat Transfer Characteristics of all Brass Cores. 
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Figure 22 * Combined Flow Friction Characteristics of all Brass Cores. 
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Figure 23. Sample Recording Trace of Downstream Temperature Response and Corresponding Derivative 




N 



Figure ?d. Comparison of Brush and Moseley Recorders. 
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SUMMAffl op HiSAT TRANSFER AND PLOW FRICTION R33ULTS 
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table III 

SUMMARY OP KELT TRANSFER AND SLOW FRICTION RESULTS 
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TAELS IV 

SUMMARY OP HEAT TRANSFER AMD PLOW FRICTION RESULTS 
STAINLESS STEEL REFERENCE MATRIX 
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SUViLARY OP HRAT i'RANSPiffl AND PLOW FRICTION RtJSULTS 
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SIMiAKY OP HEAT TRANSFER AND PLOff FRICTION RE£ iULTS 
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TABLE XII 



SUMMARY OF GEOMETRICAL AND PHYSICAL PROPERTIES 
BRASS L/D TEST CORES 

n 

1. Constants Common to all Cores: 



Matrix Material 
Specific Heat, c g 
Thermal Conductivity, k g 
Material Thickness 
Hydraulic Diameter y D H 

2. Brass 0.5" 

Flow Length, L 

Frontal Area, A,. 

Volume, V 

Free Flow Area, A 
’ c 

Conduction Area, A 
’ s 

Heat Transfer Area, A 
Matrix Density, P 

m 

Porosity, p 
Compactness, (3 
Weight, W s 

3. Brass 0.75" 

Flow Length, L 
Frontal Area, k c 
Volume, V 

Free Flow Area, A 

c 

Conduction Area, A 
’ s 

Heat Transfer Area, A 



70-30 Brass 



(Btu/lbm deg F) 


0.092 


(Btu/hr ft deg F) 


57.0 


(inches) 


0.001 


(feet) 


0.0013998 


(feet) 


0.04167 


(sq ft) 


0.05875 


(cu ft) 


0.002447 


(sq ft) 


0.051698 


(sq ft) 


0.007052 


(sq ft) 


6.9650 


(lbm/cu ft) 


61.7743 




0.87996 


(sq ft/cu ft) 


2846.34 


(lbm) 


0.15302 


(feet) 


0.0625 


(sq ft) 


0.05738 


(cu ft) 


0.003586 


(sq ft) 


0.050537 


(sq ft) 


0.006842 


(sq ft) 


10.2037 
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TABLE XII (Cont) 



Matrix Density, P m 
Porosity, p 
Compactness , J ^ 

Weight, W g 

4. Brass 1.0" 

Flow Length, L 

Frontal Area, A,, 
fr 

Volume, V 

Free Flow Area, A 
’ c 

Conduction Area, A 

s 

Heat Transfer Area, A 
Matrix Density, P m 
Porosity, p 
Compact. less, jS 
Weight, W g 

5. Brass 1.5" See Figure 10. 

6. Brass 2.25" 

Flow Length, L 

Frontal Area, A,. 

fr 

Volume, V 

Free Flow Area, A 

c 

Conduction Area, A 

s 

Heat Transfer Area, A 
Matrix Density, Pm 
Porosity, p 
Compactness,^ 

Weight, W 

s 



(lbm/cu ft) 


57.936 




0.8807 


(sq ft/cu ft) 


2845.43 


(lbra) 


0.20776 


(feet) 


0.08333 


(sq ft) 


0.060833 


(cu ft) 


0.005069 


(sq ft) 


0.05378 


(sq ft) 


0.007052 


(sq ft) 


13.93 


(lbm/cu ft) 


59.2523 




0.884076 


(sq ft/cu ft) 


2748.08 


( lbm) 


0.30035 



(feet) 


0.18750 


(sq ft) 


0.060625 


(cu ft) 


0.011367 


(sq ft) 


0.05357 


(sq ft) 


0.007052 


(sq ft) 


31.3425 


(lbm/cu ft) 


64.8924 




0.883628 


(sq ft/cu ft) 


2757.32 


(lbm) 


0.7376 
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TABLE XII (Contd) 



7. Brass 3.0" 



Flow Length, L 


(feet) 


0.250 


Frontal Area, A, 
fr 


(sq ft) 


0.060625 


Volume, V 


(cu ft) 


0.0151516 


Free Flow Area, A. 

’ fr 


(sq ft) 


0.053573 


Conduction Area, A 
’ s 


(sq ft) 


0.007052 


Heat Transfer Area, A 


(sq ft) 


41.79 


Matrix Density, ft m 


(lbm/cu ft) 


60.5958 


Porosity, p 




0.883628 


Compactness, 


(sq ft/cu ft) 


2757.32 


Weight, W 


(lbm) 


0.9184 
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TABLE XIII 



SUMMARY OF SLOPES OF BRASS TEST CORE f AND j DATA 
WITH CORRESPONDING L/D„ 



L(inches) 


l/d h 


Slope (f) 


Slope 


0.5 


29.8 


-0.885 


-0.45 


0.75 


44.7 


-0.915 


-0.556 


1.0 


59.6 


-0.915 


-0.607 


1.5 


89.2 


-0.965 


-0.832 


2.25 


134.1 


-0.941 


-0.802 


3.0 


179.0 


-0.900 


-0.796 
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APPENDIX A 



Description of Equipment 

The equipment used for the transient testing of core samples con- 
sists of a flow straightening device and air heating system which pre- 
cedes the matrix test section. Pressure taps measure both static pres- 
sure and pressure drop for the test section and the orifice for flow 
measurement. Thermocouples measure the temperature response of the 
matrix exit and the air temperature at the orifice. An ASME standard 
orifice is used to measure the flow rate, and a prime mover provides 
the air supply. 

Air Supply 

Air is used as the working fluid and is drawn through the equipment 
by a 30HP, multistage, Spencer Turbo-Compressor, which is rated at 550 
cfm using 220 V a.c. 

Flow Metering System 

An ASME standard orifice section using D, D/2 taps and thin edged 
concentric orifices was used for flow measurement. Orifice diameters of 
0.775 in., 1.232 in., and 1.540 in. were used with a 3.08 in. diameter 
pipe. 

Air Heater System 

The air heaters are bakelite frames wound with .0031 in. diameter 
nichrome wire spaced 1/32 in. apart, 50 to 52 wires per heater. There 
are two heaters per frame and fourteen frames totaling 28 heaters, which 
are wired in parallel across the input voltage source. Each pair of 
heaters is controlled by an individual selector switch so that depending 
on the flow rate, sufficient heaters may be used to give a twenty degree 
temperature rise to the incoming air. 
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Matrix Holder and Test Section 



The matrix holder and test section are constructed of closely 
machined polyethylene plastic; the tight fit minimizing air leakage. 

This section consists of a casing and a sliding drawer to hold the matrix 
under test. The casing has a removable frame on which the upstream thermo- 
couples, are mounted and the upstream and downstream static pressure 

taps. The sliding drawer contains the matrix and a plastic movable frame 
with the downstream thermocouples, t^» The flow channel is 3-1/16" by 
3-1/16" and can hold matrices of up to three inches in length. Matrices 
were placed in the holder and were surrounded on all four sides by styro- 
foam plastic insulation. This insured a snug fit of the matrices and 
also lessened heat loss from the matrix to the holder. 

Inlet Cone and Flow Straightener 

This section was designed by Piersall (13) and provided a uniform 
velocity profile to the air entering the heater section. 

Pressure Measuring System 

Pressure taps are located in the matrix holder upstream and down- 
stream of the matrix and in the pipe at one diameter and one half dia- 
meters on either side of the orifice. Each pressure tap is connected 
by flexible tubing to its corresponding manometer and draft gage. The 
following instruments were used: 

1. Ellison Draft Gage Company, 0-3" inclined gage. 

2. Ellison Draft Gage Company, 10" manometer. 

3. Ellison Draft Gage Company, 20" manometer. 

4. Merriman Instrument Company, 120" manometer. 

5. Precision Thermometer and Instrument Co., mercury barometer. 
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Temperature Measuring System 



Temperatures are measured in the system at the orifice, at the 
inlet to the system, between the heaters and the matrix and at the matrix 
outlet. Thermocouples are used for all measurements. The thermocouple 
at the orifice is a single 30 gage thermocouple, referenced to an ice 
water junction (i.e., 32 deg F) , and output read on a Rubicon portable 
potentiometer. 

There are four sets of 30 gage iron-constantan thermocouples, each 
set consisting of five thermocouples connected in series. These thermo- 
piles were constructed and wired by Traister (17). Two of the four sets 
of thermocouples had each thermocouple individually wrapped with teflon 
tape to prevent shorting and were placed in an aluminum tube mounted in 
a frame at the exit of the inlet cone. The aluminum tube acted as a 
radiation shield to prevent the heaters immediately downstream from affect- 
ing the readings. Each of these sets were designated t^ and measured 
the temperature of the incoming air. The third set, designated was 
placed in a frame between the heaters and the matrix. This set was 
"bucked" against one of the t^ sets so that the output of the two sets 
measured the difference between t^ and t^. In a similar manner, the last 
thermopile was wired to the second set of t^ and measured the dif- 
ference between incoming air and the matrix outlet, t^ - For data 

taking, the output of - t^ was led to one channel of a Hewlett- 
Packard, Moseley Division, Model 7100B dual channel strip chart recorder. 
This thermocouple output was also led into an Astrodata model 886 Wideband 
Differential D.C. Amplifier where it was amplified 1000:1. 

Differentiator 

The amplified thermocouple output from the D.C. amplifier was fed 
to the differentiating circuit on the analog computer. The computer used 
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was a Donner Model 3500 portable analog computer. The circuit used is 
shown in the theory section. The differentiated thermocouple response 
was then fed into the remaining channel of the strip chart recorder where 
the derivative could be compared directly with the undifferentiated signal. 
Heater Power 

Power for the heaters was supplied from the 250 V D.C. source in 
the laboratory. For low flow rates, the voltage at the supply panel was 
adjusted to give the proper temperature rise to the air. 
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APPENDIX B 



Data Reduction Relationships 

This appendix summarizes those data reduction relationships of 
importance in calculating flow friction and heat transfer values from 
the collected data. 

Geometry 

Accurate determination of dimensions and physical constants are 
necessary to minimize error so that comparisons between cores may be 
made. Compact heat transfer surfaces use three geometric parameters 
that allow comparison to be made between matrixes. These are: 

1. Hydraulic Diameter 

4 x free flow area 



d h- 4r h 



heat transfer area ^k-^L/k 



2. Porosity 

free flow area 

P = 



= A /A, 
c fr 



frontal area 
3. Area Compactness 

heat transfer surface area 



A- 



matrix volume 
Dividing (B-2) by (B-3) results in: 



« A/(A fr L) 



(B-l) 



(B-2) 



(B-3) 



r h = (B-4) 

Through this equation, knowing any two of the above relations leads 
to the third. 

Maximum Slope 

The following sketch of the generalized cooling curve will serve to 
help explain the maximum slope relationships: 
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Given that 

N tu CX7ld A. 



where : 



the maximum slope of the above curve is a function of 




T " generalized time variable ^ hA0/(W C ) 

s s 



N - hA/ (me ) 
tu P 




and: 




Furthermore 

tf, - ti 

and its derivative is: 



W s c s 



de 



— ~ t£s + / 

tf, - ti 



(B-5) 



d 




/ 

(tf, - ti) 



d (£f 2 - tf,) 



(B-6) 



Combining equations (B-5) and (B-6) gives us: 
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= WsCj [ cl (Ji±i z tlx ) 

'WCp cL& (B-7) 

Max 

From the generalized curve: 




x/ (chart speed) = dfi 




d(t f2 _t fl* 



= Y 



t -- -t . - m 

f 1 i 

combining with: 

c /(ic , . matrix capacity c /c sec 

s s p flow stream capacity rate s f 



-1 



and equation (B-7): 



K t f , - tiJ 
d (~~~ \ 

[ tyj J 



£ "m \X/ 



chart speed 



(B-8) 



MAX 



This value of maximum slope and A are then used to enter Table I 
or Figure 1 to get the corresponding value of N^. 

Maximum slope from analog output: 

The analog computer reads dy/dx directly in millivolts per second. 
This value, along with (m) in millivolts produces maximum slope from 
the following equation: 
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rl ~ 

Uf, - t,- 



) 



cl 




W c 
s s 

me 



hA AX 



. dy/dx 



m 



(B-9) 



Flow rate 

The mass flow rate was calculated from ASME Power Test Code (1) 
as modified by Murdock (13) by the following equation: 



m = 359 K d 



2 

o 



F 

a 




(B-10) 



where 



K = 






C - orifice coefficient of discharge (12) 

» ratio of orifice diameter to pipe diameter 

d = orifice diameter in inches 
o 

F^ = orifice plate thermal expansion factor 
Y = fluid thermal expansion factor 

= specific weight of fluid flowing = P/ (RT) assuming 
perfect gas. 

P = absolute static pressure at orifice (Ibf/sq ft) 

R ■ Gas constant for air: 53.35 (ft-lb_/lb degR) 

r m 

T = absolute temperature at orifice ("deg Rankine) 

P q = pressure drop across the orifice in inches H^O 
Substituting the expressions for K and in equation (B-10) yields: 






(B-ll) 
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From (1), Fig 40A 



Y = I - (o.*H + 0.35/8*) £- 



k = 1.4 for air, ratio of c /c 

P v 



x = 



A Pp (iy> Hg abs) 
P (i*n Hgabs) 



also from (1) fig. 38 

F =1.0 
a 

p = (Pat™ - i§: 6 ) (0.49/2 ibf/ft 2 - 

P - local atmospheric pressure in inches Hg 
atm 

P = static pressure upstream of the orifice plate in inches Hg 
o 

Substituting into (B-ll) the above expressions with the necessary physical 
constants to make the equation dimensionally consistent yields: 



™ = 5~8 9.SI -jLy dJ[l-(-4| + 0.35/6 4 )] -AIL . , 



V/-/3 



( p ^-,h) v ° z 



to ♦ 45 5.7 

Reynolds Number 

Reynolds Number is defined as: 

N R = °H G/ yLt (B-13) 

where G is the mass flow velocity based on the free flow area, A 



G = m/A = m/ (pA _ ) 
c fr 



(B-14) 



substituting: 



n r = 



-m D 



i±. 






and from (B-4), r h = p/p and D R = 4r } 
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therefore : 



N 






Substituting (B-3) into (B-15) 



,3 - 

yielded. 



d 

Afy L 



N r » (4mL)/(A^U ) 



(B-15) 



(B-16) 



Fanning Friction Factor 

The following equation and the accompanying sketch describes the 
flow system under consideration: 

Kc-H -P‘) + 2. - 1) + i ^ 

Entrance plow cope Ex it (B-17) 

eppcc t Acceleration friction k ppgct 




4 P = .s. ji. f> 

R 2g c p< |> 



FLOW 



Tz & 7z zmzz27nzzzizzzz2z2mnzzz2znzzz£n 
vfyn nni m min rTnTTrrmitiii mn±u \ 

m \ ni mmm rrrrmr rt m //■/ > v y /7rm~rn\ 

i J'/tti i It ) l > f 1 1 t >m~n mir / / /yy /*/7|77i 

r rrm n iurirrinri m/ / rrrrTi / J ///// / 1 

v/\n/j/////t/TZJ7rrrtT77T///i rirrnnmp^ 






uMiiinmnm // //////// // // i/u/j /77>7a 



I 



I 






a b 



J 



Z 




Because a gas is the working fluid the changes in pressure from 1 

to a and from b to 2, in the preceding sketch, are small in comparison 

with the total pressure; consequently, u u, and u, ~ u„. Since 

a 1 b 2 
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testing is done at moderate temperatures and pressures, the perfect gas 
law is assumed valid. From (B-l), A/A^= L/r^ and v = 1//° ; substitut- 

ing into equation (B-17) and solving for f yields, for the isothermal 
case: 



(B-is) 

where fa - ^ and subscript 1 refers to upstream values and 
subscript 2 refers to the downstream values. K and K are the entrance 

C s 

and exit loss coefficients respectively and are dependent on porosity, 
geometry of flow cross-section and Reynolds No. within the core (8). 

These coefficients were calculated from the analytical expressions derived 
by Kays (8) . 



C=^ Contraction Ratio = 0.611 + ,045p + .344p 

K, = Velocity distribution factor 
a 

1. Circular tubes 
laminar flow 



5.7 



K, = 1.333 
a 

turbulent flow 



f = 0.049N, 



-.2 



(Fanning friction factor, circular tubes) 

(B-19) 

Gap, laminar and turbulent flow 



R 

K , = l.Q9068(4f) + 0.0588 V 4f + 1. 
d 



K d (gap) - 1 



= 0.6 



(circular) - 1 
3. Square, laminar and turbulent flow 



(B-20) 



K, (square) - 1 
d 

K, (circular) -1 
d 



= 1.17 



(B-21) 
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4. Triangle, laminar and turbulent flow 



K.(triangle) - 1 
<* = 1.29 

K, (circular) - 1 
d 



(B-22) 



Using the appropriate for the passage cross-section and the 
appropriate flow, 

(B-23) 



K = 1 - 2K,p + p 
e d 



K = 
c 



1 - 2C + C (2K, - 1) 
c c d 



(B-24) 



Using an order of magnitude approximation, for small pressure differ- 
entials the first term in (B-18) is the greatest contributor to the 
friction factor. The approximation: 



**1 + **2 _ P P, £55 P„ reduces (B-18) to: 

2 — m i i 

£ ' [Hc/°^ -(K c+ Ke)-^0*P a )]^ 



(B-25) 



By substituting (B-l) and (B-14) into (B-25) it can be shown that: 

f p/fi 

Colburn j -Factor 

The Colburn j -factor is defined as: 

j = N N = 2/3 (B-26) 

J St Pr Gc Pr } 

P 

Substituting (B-14) for G and multiplying by A/A yields: 



Ac 



HA . ♦ N b 2/3 ; but N = hA/ (me ) 

A Pr tu P 



therefore : 



J « N. . (A /A)N 

tu c Pr 



2/3 
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and by substituting equations (B-l) and (B-4) it can be shown that 

j oC v/fi 

Heat Transfer Power and Friction Power 

Relative performance of matrices under comparison can be determined 
by an evaluation of the heat transfer power vs. flow friction power. 

The higher the plot of h gfcf j vs. E ^ the better the core (9). 

The heat transfer power per unit area per degree temperature differ- 



ence is : 



K= ^'3 

Evaluating c p ,^& , and Np^ at standard conditions of 

500 deg F and one atmosphere for convection (9), 



c p « 0.2477 Btu/lb deg F 

JU - 0.0678 lb/hr ft 

= 0.0413 lb/ft 3 

N = 0.671 
Pr 

Equation (B-27) evaluated at standard conditions becomes: 



(B-27) 



h . . = 0.02195 
std 



(A )( A/ «V 



3rd 



hr 



(B-28) 



The flow friction power per unit area is (9) : 



(B-29) 



Evaluating equation (B-29) at the Standard conditions shown above gives 
the equation for 

7 [aj f < b - 3 °> 



E = 1.11 x 10 
std 
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APPENDIX C 



Digital Computer Program for Data Reduction 

The digital computer program used by Traister (16) was written for 
a CDC model 1604 computer. The acquisition of an IBM 360 Systetr and 
attendant removal of the CDC 1604 computer required the program to be 
converted from Fortran 60 language to Fortran IV language. This program 
takes the sample core geometry and all raw data and calculates all the 
heat transfer and friction results used in this report. 

The program uses a curve fitting interpolation subroutine to deter- 
mine the value of N from maximum slope and conduction parameter. This 
subroutine uses Howard's conduction parameter data. (Table I). While 
not used in this investigation, a cyclic method for finding the value 
of N tu exists and a subroutine for this method is included in the program. 



95 



LU 







— 1 


CO LU 










X 




UJ 


< 


UJ o QC 










•> 




o a 


•o 


O K — IOOO < 










UJO 




2 — 


CAK — — 


< <— • 










CL 2 




O —COO 


— LU <00 


CA QC X ^ O 










o< 




2 KXZ 


LUCA — >~ 


LA < 2 QC LA 2 










-1 




< U-O— 


< LL — < X 


< — lOlL II < 










oo> 




— LUO 


K CAUJCL 


CL0C3 — XZ QC 










#> 




Z KQO 


<2co— QC 


<OKKXLA UJ ► 










XX O 


OlA 


OlA 002 


O— <UJ 


X-IOC< 3 — Q.Z 










3 • K 


LA * 


— o <uj<o 


2 LA X 


ODhXXJ 3 










XLLZ 


— O 


KO IOKK 


oox X OCK 


— JOOQC — OOooqC 






— 




— QC<OA 


LL *— I 


ox — < 


— K32LU 


LL 2 OCCOZZ 






o 




X<CLK 


»w 


-hOI ox 


IOJ2U.00 


< < U-K <DI 






UJ 




< 002 


— — 1 


OC UJ UJ O UJ < 


K2O3002 


UJ— 2 < • QCO 






CA 




X IK UJ 


OCA 


U.2K- -JUJQC 


UJO— 12 — 


XOCQC— X — -1 < 






X 




O X 


LA ► 


LA IK OC O 


• -J o<< 


• KKO X<-1UJ 


• 




00 




UJ— CALL 




DcotuDr- oo 


O *OQC K Z ILUJOOO — < 


QC 




K 




ZlXh 




O — K K O QC 


uj3:k K 2 X LU < x— — OCKOC 


LU 




O 


— * 


X — 3CO< 


CC LA 


-J-J 00 CL 


K OLL • 030 OK UJZC 


K 




o 


CA XOX -IK 


»w> 


lloluqcocll 


LAO — KO.J QC • LUK K LULL 


UJ 




> 


CLKZKQC ►— JCA 


— 2 


>-0DUJOOUJ 


LU U_ C?QC < 1 


o>a>txiu<<x 


X 




— 


OtA< OUUZO 


03 


OO K X 


h- oOOOJQOQtU. O0CliJ2:w< 


< 




— 1 


— UJO 




IAK 


2 OO-JOCK 


oo clXqc 


t- ootxx art- 


— 




-J 


KOOIIUJCL 3 


w ► 


— CK <UJLU 


UJ — 2— <OLUOZlLK OUJ< 


o 




w 


o 


Ooo2 Cl 


LA — 


22 OCCOCG Z 


qC — UJO— >2! •< OUJ 0.0 






X 


OO00O K LU 2 


► vO 


— <UJ<<X — 


C— '- >*o 1 


LUO— 02 — K X 


UJ 


• 


w 


X — 2 


• O ►! — 


— — 


X ctIJDl 


O 00>-«3 Cl'-'UJ CL LU 31 IX. < LULU 


a 


— * 




t-jroi- h- 


O — 


QC300 2K 


Z< — LXZLLO II «>-Z> — C£ X 


— 


CA 


K 


uliod 


a. lu 


LA< 


UJOO UJ — 


ujx— coa:0 


OSJOOILh 


a 


K — 


3 




4-rax 


—a 


K— 1 Ol>3 


ID 0300UJK aOUJU-O 




—1C 


— a 


OOZOUJXK 


UJO 


LU CD X — h- 2 


K-J<QCOOUJ3rXUJ<XOCLQC 2 


o 


O AJ 


— O K 


o o 




o ► 


C OQC <- 


OujO ookkoo O caujqC— 


K 


• >x 


O AJ — 3 


— OC 


Cl 0_ O < O 


o— 


UJ<LL)00 


QCOQCQ.CA 


a:o • cllu< 






AJXooO 


-IOC000 K 



*o ocoujx— uj n 
— IA OO 3 



X k ^ou-o:hh * • 

ZJOO 2UJUJUJ 



O-J 
LAO 
— O 



— — O 



O •< •—LU ►COLAQ XX COK LUX JZ 
LU CA l>0 XUJUJ — U.UJ2- 0CX2 OOJh 

O — U.ZQCZCOJ • ZZm>o 2c 0 0i<0iODO — 2 — ZhCOIOOO II UU2TOC 
LA OC — 0 I £*2- 0<X0 O^^CXJhU-ZO LL — X — — >— I —2 O*— << 

— ► O00LL<<< UuKO — ^Z OCX • C0< — — — LU — w<ZZh<OZhQCZ 

O — -J U.OlOk<jh 2 •OoLUf^oo XOUjZcaqCoC O K-JZ*— 0< »<Ch 
UO COLUO_| CCOOOCOZXrZ^lU ho LUX UJOUJUJ — — XLU-J< — — QCOC-hK LU 
►LA 3IUJ<0A2<UJ OUJX— l<E<LU — ►ZDtLUCAKII CL OCLL □ — «J — K LL XKcD 
— — Zha M<oh>- 000—1 OoCUUUJKX3 ooil ID — AJqCZXX *<CL ►ujkqc 

<11.10 OC o JO o UUO»— K LA3ZO00K0CXK — — OLA— JXUJCAX — <2 
zzaoi- UJOhJOOO ZZhffi<OZ O — C — Q<C < QC oo_JUU>OK3:I UJ 

<ouj xiziuujcoo •— x<xoo — xoo— ac — ozac zxxzu.ujh k< ox 

O ► OCQCClCK *<K — X — O K DID^U-Dm h^I — OLUCA — — — qCO UJ|IIO K 

— O LL CL LU h- Cr: UJOULOC • LL — -JK 2<LL_J < K QC K O CL CO — LA CC UJa^hJ KUJIJJ 

o oo<ozcouj o k 20 CLuua^cr:^ uj<ujxo ooh-oa olu o— oik 
K OlA CL 2 LA xoo 02 y-OJoOW OO X >t K QC 00 UJ oC *0<Uj0o -JUj #:x<KUO 
1 law clo— < oroc oo uuo la— qc — ujo zoC3 lu — koooclXoc-i— cl 23— _j o 
O — O — OO.QCUJUJK — 2 00O<IaC UJXUJ— ICAX ► <CAO OOAqCX— C? OCL2X 

► Oh OK LA CAUL LU QC X J*— 2 LA —JO K < LA >0 K O O <C0 O ULJ LUCK LA UUUU — LUZLUC) 

X ► ► Z<C — ojXOK ►O — <LU 3U0<_JUJOQC — OUJOCLOC CCXXK02LA*— QCCtCA 
I — — — XaC ^120 00 — X*— O *OOUJC2 <0 •LU'— an UJOa K X OK LL LU 

<200 LA0C<O< O LA 2 X— LL2XLAJJOC «0 rOQ-LLC CL XX 00 LL O OOLJJ O 
wOLAlA OOOQCO ZUL2UJ00QCOX •— < QC — J UJ^oO — KUJ< — O OKUUOUJ C£LA2 

OOK— w — LL <AJ • C — > — — < — 2>aC<K • •<< OO QC-JOKO — XX UUK2 K<Oj< 
* ►ZXOaCZCO O <UJ LA — QCLUUUO *— — a-*— ••— OKOUJmjhKkOj — XX 

—1—10 ►Aj O *— K ZUJLUOCKlACMO KKXOQCXUJX QC- LL QC OILUJK — XCLUK3CK0 K hUO 

CLAQC LL <<— KILUZ — << 0COUU2LA2UJ2— || <K — OC QCOO^QhZ 

uu ~XO_J o:OK LAKC0OO2OLUX ►CCLLL o — OK — QC X JOOC KZXaLLOJUJJ ►— O 

QCC LA UJ LA LU < ►LU XOZXCQC UU O > — < OCOOQlO XlACOX CQKCOO 2 

ClZwCQZLLU-OOKOO 00qC<<A>lA2022KKX0 K >— I —I CL >* K QC QD — Z — 

KOOUJCOCAO 00 ZZAJX— 1 #<X*— O LUOLUOOujOLUCOLACjCLLQCLUOUJ ►•— O 

— — --I — 2 K UJ— COO;* LL <AJ 0>OOKOOOKOCLOXK -J — — ICX 2 

0200— ►K < X CALL QC OLU • O — OK — X ►UJ_| — O— CA II II If II O-JO OC 

— 020— O0CO0COO2XOO —I— I 2 UJOOKZJLU II II || II II II II II II 202 00 UJ CL 

OXUJlA — OKUJ— 0<KZUJ *OLL2 — ZLU020000K O^AJCA — > — 20C 00 

O-XX— LAqC OLLO X O — < X 2I« 30CA — » CM (A >cMA O K CO O'HH^^oOOIA X OUJ 

XO — LU — KK JIOIO- -JqCUJOCOKOqCXO LU OLAqC 

— OOoCKLACUJUJUJO KqC>XOLU— JKJ OAK O ZUJ qCQC — I _J — J — J— J — J— J — I _J -J— J — J — J QC QC OC — J < QC 

•—AJ2UJ2I — <>* — < — XVqCOCOOCZqCOIOOOOOOOOOOOOOOOOOOOUJO 
— I0KU-UJCD3O0- lU.OOOKQ.UJU.OKU.OOOOODUOUOOOUOlUUJ.o2o 



oooooooooooaooooooooooooooooooooooooooouooo 



96 



• • 

o 

<■ 

LU 



a 








CO 






















CL 














o 








► 














l- 








oc 




















z 


LL 


► 












o*» 






► 


► 


X 


► 










LUO 






o 


CL 


a 


X 










o< 






»— ■ 


X 


>• 


H 










ZUJ 






K 


K 


o 


► 










<QL 






< 


< 


* 


m 










X 






ac 


► 


z 


• 










uo 






► 


X 


< 


v0 










h- 






X 


o 


CL 


LL 










UJ 






o 


> 


c0 












CD O 






► 


c 


h- 


x<to 










LU 








► 


* 


•— < • 










oo 






00 




00 












00 2T 






► 


< 


X 


coll 










o< 






X 


a 


► 


• ► 










XX 






o 


00 


X 


SOX 










o 






► 


H* 


-J 


LL *— < 










cr 






X 


► 


UJ 


► ► 










O 






UJ 


00 


a 


xcvi 










h- (\J 






LU 


X 


► 


• 










z>o 






► 


► 


o 


►ip 










UJ 






X 


X 


a 


P01L 










ri- 






3r 


— i 


► 


• ► 










luz 






— » ► 


UJ 


c 


IPX 










L-UJ 






^< 


O 


a 


LL •— < 










<x 






s ► 


► 


-i 


► ► 




o 






K— LU 






o 


o 


UJ 


X CSJO 










UOL- — 






< 


CL 


o 


^ • • 




4* - 






< — • 






► 


► 


► 


►IPO 




CVJCO^lPfOh- 






•4- IP 


X 


2" 


CL 


o 


K 


OsJU^vO 










Z /> ► 


3: 


► 


a 


a 


-J 


• ►* 




► ►IP ► ► ► 






c -» 


► 


O 


CL 


-J 


o 


IPXsO 




r-opj-tmosj 






mO • II 


< 


M 


► ►► 


LU 


> 


LL^CP 




PJfM ►PJOsiCVJ 






L-Z — ~ 


► 


H* 


< XOsj 


O 


o 


► ►CO 




► ►IP ►^ ► 






w<t/0 - 


o 


< 


h-CL — 


► 


h- 


Xr-H • 


o 








o o — 


< 


CL 


< ► ► 


H 


► 


-• *o 


Xr-. 


OsJPJ CMOOsJ 






o * o 


► 


► 


Ooo m 


-J 


00 


►IP\ C\J 


xo 


*— .^«*"*«**o*»- 






<— — r-4 


cc 


X 


< • 


O 


O 


IP ll o<^*— < 


LUO 








zr - 


o 


o 


h- -o 


> 


► 


• •* ♦OO x 




oooooo 






2T<wvO^— i 


CL 


► 


3-I-* 


o 


< 


sOXlP •< 


* ► 


CD OOO ooo 






II 


► 




CLXUL 


h- 


L- 


LL O • 


<o 


1-0000 «o 






O0w ► ►-* 


X 


00 


ZXiP 




—*LU 


•> ►* OO 


1 o- 


* OOOO 1 o 








CL 


► 


W ►S 


OOO CO 


xc\j— a* 


OO 


toooo-o 






00* — II II — 


► 




X-|vO 


o 


• ► 


^ • CO 


r- r-^ 


1 • • • *f\l • 






oJIUOmh^ 


00 


o 


— iUJ • 


►mo 


-IPO* — 


NO - 








ou ► ►- 


-J < 


►««■* 


•— 4 CO O 


<LLO 


COLL vOZsO 


^O 




O 




Z*-*00 — ^--5 


LU ► 




► < f—4 


H • 


—4 ► 


• ► • < • 


'to 


sDlP vt OOlPIP 


o 




<h-* W~w 


CQ -J 


UJ*— • 


X JIL 


UJ •— < X 


iPxocLm 


OO 


f*H so ip Nt pj ^ 


o 




XOOOww^J— . 


< X — 


LU -■ 


— O 0- 


CO 


’ ► 


LL ►H 1 COi— <h“ • 


>t • • • • • • 


OO 




O wQ_ < Z CO O 


-J xo 


O't't—V. 


►o— 


► — h-h-^ 


• O o 


'O | | ll II 


OO 




v.<az)w • 


— • 


— * • 




• 0- 


xmjwOo>4.^ « 


•<<<<<< 


OsOO 


o 


LUCLOOV- win 


— >t — OOoO^^< 


►oo^ 


— * *0\CL IP , ~ m4 O 


OLkHh-^L* 


O^sO 


o 


OOXwww ll 


h<o^O^^IOO h i— <ocd 


i -'O > X w 




* UULU LULU LUUJ 


IP 0 s • 




LUUJ — — nO 


O O O LL <f LL 


i| — u.-^cx:fOij.oQ l 


-t-cro + o 


IP CD CO CD CD CD CD 


OsJ^O 


r^- 


Xh- l| •— * »— ■ < 


-nosjcnh- 


►in O'*' ►*■■» ' 


xr-fM ► ►•— ►>— vo-om *ooocr www^^wopjo 


00 — < 


o C - - *w 


►* — - 


in w 


w sO 


►w>ca: 


wXwG3!CLOOa^(lC0 00c0 00 0000l000C> • II CNja 


— LPinirh- 


in h~ ip 


wh- 


h^i-tin>-^u.kH n hi • i x 


XOCDXICOCCoC O 


IPO 


IP 



ii w<ao w<^< <UJ<LU<OW<UJ ►< ►o n <luoo:ujouj *<<<<<<o • n xo • 

o wOOOX OXOXOXL-Xl-Xr-C>Xl-CLXfnCL wf- II ULhkhOOOOQOO>-OX3hO 

LUO 00 <<<<X n ^-koc<oc<oc wart-no: *so n 11 o:^ h n h ww^www m h 

-i— I UJUJUJO H CUOUJOUJOa:OQt:OCUJOQc:H-OvCUJ_ir0OLLU.CDOcClXLLa,LLtLLLU-O<MUJ00a^ 
LLCPOarot'a'lL.-JXOCLL OClLtt'LL 3lL3:iLQCtU.?<a.U.I-00aL-D^KOI-riMwwMi-.«OXO<OX 
UJ OsJ ^ 



accr o 


«— < 


C Hi ^OOO 


>t 


IP^ON O 


p- COO 


o 


^OsJ^CVJO'tlP'OP-^ 


OsJ 






oomoNto 


f— 4 


O ^ 


f-4 


o 


O CD —4 f— 4 f4 f-4 fM f-4 OsJ 


OsJ 






ro ro <t 


o 


o>t o 


o 


o 


OO 












H <“ ( 




•—4 


F^ 





ooo 



97 



CD 

►— 

h- 

44 

CD 

44 

CL 

O 

o 

* 



o 

+ 



LL 

44 

O 



m 

CL 

44 

O 



O 

X 

o 



a.<z 

-JhaJ 

LULU^s 



X 

3 

* 



cl 

a 

00 

O Q 
O' 44 
O' < 

r-4 00 

•^00 















. — . 


< 


O CO • 


O 


< 


<<~ >^OH(M 














< 


* 


w-vo 


a 




OCMOOOO 














44 


44 


b»^0 


X 


'V 


o • »ooo 














44 r- 


< 


occxro 


UJ 


— * 


CMOOCMfMCMm 














<< 


K 


OLLO 


h- 


— * 


•* | 4 • • 














KO 


UJ 


oODh 


44 


X 


<wrsjmvOr-o 














LUsO 


ao 


OX.W 


< 


QCW 


O* LLOOOCM 














com 


44 


44—4- 


1 


LLOQf 


044 44000 














l o 


m 


oxx 


o 


044 u. 


cm— Ocmcmcmo 






o 


o 






• • 


m 


OwO 


m 


SJD 


y • ► ► •►CO 






o 


o 








• 


<44 a oo 


m 




— — <msof^ * 


oo 


oo 


oo 


mo 


oo 


oo 


w 


+ 


•O 44 O 


m 


xxx 


• 4-44 OOOO 


oo 


oo 


oo 


OCMO 


oo 


oo 


h- 44 


f— ^ "tl - 


CM 


w* oooacooooooo 


*-«o 


OCMO 


or-o 


0^0 


omo 


OOO 


CL 


< 


44 >DX 


00 


0044 OO^nOcmcmcm 


oo 


OCMO 


OvOsO 


OCMO 


OO'O 


omo 


OO 


• 


II H0>J 


• 


o* •«- 


moo 44 — — * 


>00 


<>00 


0000 • 


moo 


mmo 


HflOO 


ooamo 


ao •uj 


1 <QC''tIf f >^0'^C'O'O'(7 v 


-<>0 


00 ^ \0 


<rHO 


00 »>0 


^-cm>o 


OrOvOXO JC l — o *<0 


h-^<o— — 


m i <o • • • *o 



o *00^0 *C0<O 00 O O • CO loo *C0 vDO •O^LU'v | ZO'*-* 4- <UUhoN • t— • — *0 •rOCVi^O^: 

•Ccvio *o(\j o • ii (nj lo ocnjo #ocmo #oi/) »oko ao iimo:do< *m ii o£*-»x • «-* t i i i x 
on mo ii *1111 soo ii voo ii oor-i k #^-»lo xozoo* o ii ti —o 

UIO *11X0 HIXCO IQ *l|XO • II X< II II o II U*-* *CL CL~— II I— II UUUJUJUJ II 

XDhOXDhOXDh tl XDhOXDhOXO II II ll II m-* X< IlhOd II XXX CL II UJUJLULU 

-Jt/5 II -JOO II -J00 00 100 II -J00 II -JOOaC f— CCD XOQCQCZIIII 3—^0^ o*-~ — ~0 

LU <S) O CM UJ lo o CM UJ OO o X UJ C/0 O CM UJ 00 O 00 UJ l/> o ^ I- 4- 4- O — Z O O O OC ■-« D* Z 4- IL 5*: Li. NJ X LL U. LL U- x 

ocoxocoxcco Q<oxc<oxo<ocohi-hQCa:uooaooiofQrxHiixN.MMMX 

ro < m nO r>- ao ^rgro o (\j^o m>or w f'-ooa' <MT\sor^o 

rvj CM CM CM CM CM CM rOCOsO h^hO'O'O' OOOOO 

OOO ro vOnOsOOOO' OOOOO 

CM CM CM CM CM 



98 



UJ 

a 

o 



oo 

UJ 



m 











# 






< 










* 






QUL 










— 






X\ 










• 






<0 










o 






•J 










o 
















o 
















r—4 






< 


«*» 














c 


o 








— 






X 


o 








X 






<oc 


* 








— - 






JLL 


o 








UL 


o 






o 






— 


a 


• 




ac 










w 


o 




X 


\ 






r- 




o 




or 


«*■* 






• 


i— * 


vO 




t—4 








O' — 


X 


m 




< 


• 






ir> — 




* 




z 








>to 


UL 


a: 




o< 


^ l 






+ * 


— * 


IL 




UL 


•— 


Q 




oo 


Xm 


* 






ire 


X 


— * 


CL*~ 


C£* 


rx 


— 


UL 


* X 


X 


a 


x\ 


** 


m 






* x 


► 


o 


UJ — 


• — - 


00 


X 


K 


O'* 


o 


«* 


KO 


nT — 


(XI 


<w 


c 


O • 


X 


— * 


*m* • 


— X 


• 


X 


OK 


a os? 


X 


X 


^ st 




vO 


K 


XX 


<\j 


*• 


— ■ 


cx>fr 


— * 


* 


► 




* x 


— o 






X 




a: 


* UJ 


X* 


m* 


o-t 


o 


X 


z 


vf~C£LJJ 


— '-J 


• Ql 


* w 


* 




z 


omoO 


- o 


ox 


ro-J 




crx— 


00 


vtX 



(\i «oa x 



• "vXO 

o 



— m 



itmli 

wO 



io ^ 03 : 0 - 



-j<a: 

ujo 





o 


oc^ua: — o 


O'sh 


• O 


x^ X 




K 


»— (\X 


xx— coevx 


• 


1 


1 OUC^X'O 


J — < 


com 


— xa 


-J 




vO— UJ 


wwX — < •‘O 


+ 


o 


oa^aoxro 


oouoa 


aoo 


ooo\ 


X 


* •»* 


•*st ► 


UJO^X-JXfvj 


o 


X 


x* ♦ * — U*>\ 


* <* 


# sO •rotxi 


• • X Q-CD X 


inr- 


*— 1 (XI CO 


QhLL^ OHO 


X 


X 


xmoc^o *a 


O# — 




CM>ta:o 1 


1 


vO *a: 


O* \Z> — <XJ 


X 


* 


* 'tOXO'OO 


CLXX 


a* voo 


\vT * UJfVJXcO^Q 


— * w 


• — ^-—Kin •* 



O* OKOO'O* XO~<* X<0~ CL — *-HNtOULC£*(X.~» ^H~UJ c£oX* 

r\i'0<X'-H(XCX • • # r-m-o'jwC +<IO • -J<x* *#h#c^OhXOJ:j (XI UJ — Xf^X 

O «0 *0 • + <XI • ~LnX*«- o O'-'-J •* 00UJQuXNf(XI I mo • • II 0~0 OK II JO ‘ 

| (XJvO^^X II M X •XUJOO.Q.Q4' Q.'flUHU.OOH^ • OKQ-KZ OKsO^ — 

— • • I — KO^ II X I | ULLL II ^ II X I — O II •> ^ 

O II C M O II *UL< II || — X + X II O0K<^Q- It II UL II K II -~GC — O X II II UU<\ \ 

K K K ,^HhZ^O0^X^^IKai<Q.a-'X|l II ^tOl^ZLL l| JXhJ^^^hJtS^ 

a o a ii ii ii ~q£xolX— luo a— i n n n n n o^n^izqid- «j~ jjir-c^xh- 
axoxoxoxxa:o~<~oo u cx ii uuo<^fvjxiu.Lm.^OhwQiLLroo<o^a:oxz 

0X0 XOXOOUJ 1C U. 00 (^ 00 XOOCQl a Q. a CL a: ULU.LLLLOOOUJX'-'UJOKOOOQrX^U- •> X 

^rvjm 



»— i 


(XI 


mo 


in >o oo o O f-H cx 


m 


-t 


LT\ vOK 


O' 


HINJfOxtvt IT\ 


o 


o 


OCX 


0'0'»-H-HCX(X(X 


CX 


(X 


rxfxosj 


>o 


>0OOOhh 


a 

(XI 


o 

(Xi 


oo 

(X »-H 


00'0'0'C>00 
> % 


vO 


vO 


>c »o >0 







99 





»• 






-J 










-J 


in 






o 










o 


• 


< 


oc 


o 


s 




— 




a 


O' 




z 


— 


\ 








— 


LL 


< 




— k 


in 




l\l>t 




a 


L 


Q 




** 


• 




• • 




•— • 


rg 






«w 


o 




oo 






« 


OC 




h- 


H 








t> 


o 


LU 




a: 


LL 




LLUL 




i- 


UL 


UL 




— 










— 


— 


OO 


3 


— 










— * 


in 


2” 


h~ 




O 








*— » 


• 


< 


Z 


«— * 


HH \\\\\ 


w\ — . 






■w 


00 


OC 




3 — 


h- SSSSS 








oo 


LL 


h- 




»- in 


< in in in in in in in in in 






— 


» 






— • 


oc • • • • • 


1 • • • • 






*- 


rg 


K 


LU 


— r- 


•ecoooooocoo 








• 


<XCL 


— LL 


oil h— 


IrH — ^4 


OOO0 




— 


o 


LL1 < O 


w - 


< < U. U. U. U. U. LL UL U- LL 






LU 




is j 


oOOg 


► 








Q 


LL 




00 


- • 


X 








a 


— 






— O' 


oc 


— 






— 


in 






•—LL 


— 


X 






«•— 


• 






w — 


< 


OC 




— 




O' 




X 


ourvj 


— 


N 




rg 


— 


LL 




o< 


Q • 


oc 


-*uo 


»- 


• — • 


Q 


— 




Za: 


or- 


O m 


<0C 


l-K 


— 


O 


fO 


o< 


-LL 


Cl cC 


1— — 


Inm4 




— < 


<— • 


OOu 


— - 


•“ LLoooO 


1 X LU 


>o 






— o 






—in 


oo _j\<<cl<q:ccc 


— < 


z 


— H 








— • 


< CO 




l-CL 


3 




X LL 






ULC0 


► -J 


1- 


o< 


OC 


a 


•» — 






•“LL 


o 


< 


30 






— in 






— » — 


<< 


OC 


O 


oc 




>— i • 






—in 


-LU 




2h 


o 


— ■ 


— O 






w • 


Ooc 


oc 


o< 


OC 


«i 


LJJ»— < 






LL rn 


ac< 


— UJ 


OUJ 


oc 


o 


-UL 


< 




CY 


- < 


< Xh- 


X 


UJ 


— , 


«-* — 


►- 




•*UL 


-J_J LU 


LU >UJ 


o 




— »i 


— in 


< 


LL 


— — 


X< OC 


ocoo^x 


<x 


< 


— ■ 


W • 


o 




—in 


Xh < 


<3<<XX— 


h- 




QC GO 






«•*» • 


“Z >< 


— 0CO0C0C 


LU 


z 


-LL 


z 




zr- 


— 0CO HUJ-J 


OiQ < 


►UJI- 


CO 


3 


—m 


o 




3LL 


•juloc — oc<r 


luool^:x< 




0C' 


1— 1 — 


— * 




a •* 


— X<LLXOO<I- 


UL 0C h" < 


LOI-X 


O X 


21 


wog 


h- 


oc 


xrg 


-JX >-z z 


00 — Z 




mxr- 


- 


LL • 


o 


z 


• 


X*~X0LU3:0 


zoooo 


-XX 


— f— 4 




OCO' — 


M 




z— cr 


X -tNZOOCC><«Z«' 




O— •“ 


00 


-LL— in 


•*ac 




•OLL 


-»■ QCOroLU _jiLH-aC-JUJi-Ocr>m 


ff> r—4 —4 


11 -<■ 


— — in *— <ll 




— <rg -OCOCUL >C ► JX UL — 


)-30o^ •“ •“ 


•►ao X 



m • ^ x 

— -W 

noil— z 

N“— -►— LL— '►— 3 

^LU— < -l±J< OC 

h-H^ZUM-r 

— OC •MQCIII 

co:hO®aoooo 

Q3a:uui IXlLvOvOsO — 
»-^ *— 4 t-H rvj m 

r- oo >^in 



r- f- O' O' 



ii OcoOllc Coo < 3 ^NN 

— •>— a<w>0N fc itt:uu*-iO>-Qf<Qv0v fc v^ 

O^S^S^ChUJJQ£:<QUJZ^^ 

(>^k'j\< hj<(^ooaj>a:o H 

^ujccx lu<ulxlloo:lxx<:ouj< 

CM-X l| 3: t| f-X L-X 

oacoiLOO<^or^r-r-r-r-r--h-r-r-aco 

03:uL<o^<3tULmmfOr^rororomro3a- 

^f\jrONtuONOr^occ> -h 

O'OO^cgm^ in o r— coo— *cgcno 

-IMOOOOO OO oocooooocorvj 

C S 0 V '0'0'0'C'0 O O s0 'O 



o *— <in 
rvjso*^ *m 
org 

*: h- ii 

— » I uj< o 

^ — oc 

ii llocOlucdi-z 

3CLLGOOOOLU 



100 



CM 

oo 

oo 

a:oc 

UJLU 

b- 

22 



>tLr\'Of^coao^<Mcn^uo sor-aocr o^cMrosj-insOr- 

O O O O O O <M.— I r-H rH fH .-4 *-4 —4 <\J CM CM CM CM <\J CVJ CM 

oooooooooooooooooooooooo 

QCCCQCCCQCQCCCQCQCQCOCcCCCCCQCQCoCQCQCQCcC.QCCCCC 
LULULUUJ LULU LULU UJLU LULU LULU UJLU LULU LULU LULU UJLU 



ZZ2ZZZZZZZ2ZZZZZZZZZZ2ZZ 



X 

aO 

cc 

# 

1X1 

# 

r 

LU 

# 

O 

CM 



in 



VC 

# 



in 





0 




or 




* 




2* 




Uj 






— 


5“ 


CL — 


X 


•^sj 


* 


21 | 


<— 


ILO 


• 


•> 


nT 


COX 


1 


X 1 


— O 




Uj-— 


X— 1 



— — 

O# Oh 
>X 

* 31 

OJ^OO 
XX 30 

2 3 I 

"•hJO 
X • O • 

03 

_>*-< ii — 

xx 3 
ta x 11 1—0 
3XX OJZ 



sO 




b- 






rH 




•* 






w 




Cl 






ox 




O 






1—00 




► 




O 


a •* ► 




vO 




1— 


ox— 




r-H 




9 - 


h-00r-l 


r — « 0 


•» 




X 


► ►in 


OOO 


A 


CMC- 


00 


20.*- 


0—*(M 


M 


OO 


•» 


00 OZ 


r-H 




.-ter 


0 


- -O CM 


in 


r~4 


► ► 


2 


CL 1— h~ *■ 


i^h 


*— « 


CM *-h 


•— 1 


0 ► m 


r-H »— 4 ► ► 


X 


O (M 


•> 


~X — 2 •* 


0 *-sOU0 


00 


•— 1 * 


2 


XOOvOOCM 


CM h- r-H r-H 


•» 


► fH 


O 


UJ *THh 


► *-H 


< 


i^cm 


h- 


1-2 — — 


<— H — * <— * 


CL 


O' 


► 


ZO<JC0 


O — — — 


OO 


a — 


h- 


*— < ►— CL 00 r— < 


cm— mm 


Zw 


M. 




•*0 •**-* 


► QN *— 4 —h 


t ' "« *H 


— O 


»H 


UJ< <mtt> 


0 


•*2 


— 2 


2 


za.za.-OH 


o<<< 


*-H •— 4 


^H k-H 


— * 


— *300 1 — » 


LU CM CL CL fL 


II XX 




X 



h- — l~ mo 000 r- I s - ^a.Ohh cl 

I- 1 1 0 2— • cm o t^z I | 02 vO 

o^zazn«Q.a.ia. • xmoM-omo —* 2 : 2 : *x o 

iXJLUZwhhOUOU^^ II h- II h- II l| X 

ocxx— ^2:— — ~~ 11 1-3 o 3 X 2 — ~ ti 1—300 

Dx-au.zaa.iLu.u.DLjzczozc<CLU. oujkz 
00tfro Mw L)^ w «^hDf't-iO'-'0)‘HOOU M '- ( h^Z( 3 IX 





CM<m hoho r- vO inr^ 


O 


CMC- 


sf 


OOOO *-h 1 i—4 fH 


m 


ocr 


sO 


C\J *— 1 *— 4 CM 




—• Ch 



101 



ao & o — 1 eg ro st m o r- go o o ^ rvj m in >o f*- oo o ^ cm m m o 
<Nj<Nrornf^rOrorornr^roro>t st >t >t ''t ^t^^vnvnininintn 

ooooooooooooooooooooooooooooo 

xxxxxxxxxxxxxxo'xxxxxxxarixaiflexacx 

■ jLiimuuiumjJU j mix J iJJuumujujujUMj u uniMjjujuuujijmiJUJUiiuti < 




o 



o — 

CL— . 

x~x 

ujx * 

i_ 

Zen— 

*^ w O 

0 + 
or ♦ c\j 
qc\j* 
u.* * 

* — 
ID — i — » 
Oh^ 
ZZ — 
<hX 

xx 1 

1 



— — 






LL~Z 


H-rsJ 








2 I 






IX 


~0 






1— X*— 


>w 






3~# 


*> 


— ■ 




Ox — 


H* CC — 


H 




— X 


2* O 


2 






M JO 


•— • 




X — CM 


x<co 


> 




(MX^ 


►Lb ► 


•» 




CTi 




H- 




•^rHw 




2 




O'w* 


>*-0 


»— • 


rH 


•O — 


•"X — 


XCNJ 


l.f\ **H 


CO- ■ 


X ! —O 


•* *> 


► •> 


-H* X 


W<S - 


SI-* 


>fr .*N 


UJ-.W 


<-H w 


•» ► 


•* r-«4 


KX 


2 >- •* 


>0 


sO •» 


II 2 1 


t-^co •*> 


*N* 


— \T\ C<\ 


~h- 


-J* — •* 


X 


*— * •* r— » 


h- H~ X 2 


Cl — 1 SIX 


— > 


*-— < *-H .—»* 


22 | 1 - 


(✓:<«— w 


Z — 


+ N- — 


mm-«X 


LU X 3 


2 *— < 


X— •* X 


XXr-l^ 


LUX O 


MW 


-***■’ f-H '“■"H «••«« 


*0+ + 


2 2 — 


JX 


— x + x — 


HlXh 




3- 1 


-* I x 1 x 


OCO^Z 


KmhQ. 


00 K- 


'»*■ h-*-- — pw^«»w 




Oijoooo 


2 


> 2 Z >2 X 2>2 


h-«~>-Z 


3*-2 


>•*— c 


II X— II X— < I h3h II ^HQh< It 11 x 


X— JLL'-J 


XX hh-DXLD + 2 : 1 hXKD 1 


KZXkhD 




21 2 ^X — 2 — ZHX 


«Qi2Z^Q 


r> sr—< 11 20. n 0 


UJ U- *— < LLl II CL It QCL^UJ II O^O^^lUZ 


^^GO^LU^^OV^^V^Y‘-yO^>-^^OClLi>->-fY'LU 




fOO 


—4 c\jsf cn — 1 sooj rn 






h- 


p*- 1— < 


0 



102 



r- oo O' o — < r\J m 4- ir\ o f*- cc O' o — < <n ro 4* m o r*- ao O' o — * <\i m 4- in >o r*- oo o o *-• c\» rn 4^ it 
ir> to in so o 'O >0 so >o o *o >o 'O h- h- h- r- o- h- r- r- oo oo oo oo oo oo ao ao oo oo 0 s O' a* 0 s & cr 
coooooooooooooooooooooooooooooooooooooo 
(X<ycyocaca:<xa:acc£<ya£occya£ac<yac<yocc)ccytt(xcy<*ai:cyo£acacoca:a£aa:tt:ccac 

LULUUJUJUJUJUJUJUJUJLL UJUJIDlL LU UJ UJ UJ UJ UJ LUUJUUUJUJ U-UJUJUJ LULU LULU OJUJ Hi UJLU 



O 

O 

ro 



ro 

► 

o 

o 

m 



o 

o 

ro 

LL‘ 



a 

o 



ro 














w 






ro 








Q> 






•> 








•» 






r— < 








— o 






1 

y: 








rvlO 






W 








I r> 




<N 


<r 








Q w 




► 










*W X 




1—. 








— * 


► 




w 


r—4 




. — . 


y: 


CC 




< 


1 




(< 




* o 






y: 




— 


a 


— -JO 




— 


■W 




CnJ 




aX co 






CL 




* 


yc — 


»-i i -' •> 




•w 


1 


CO— 1 


. — . 


w 


Xcr: 4 * 


— 


— CL 


— * — . 


- 1 


— m 


Cl + 


► ► >w 




-o | 


rH 


-it: 


*'■■"4 ► 


•K- ^ 


>~o 


— 




^ 1 


1 — 


+ — * 


— .•W' 


►X * 


O 


— * — 4 *^ 


— y: 


— XCQ 




LN 


x | — 


\ 


<\j wOvJOvJ 


CL — 


ro — * 


— < 


— 1 


— <x 


— * 


— < — •. 


+ rn— 




-4 | 


^Vj — 


Qw^ 


— * 


O # CLCNi 


— # c \ 1 




# — 


! ^ 


O > 


— ^ 


— * V*— * — * 4. — — » 


»—4 — » »>— » 


| — — 


— (X) 




— .CC - 


^ — 


C\j •—* < 4 . \J 


l -h^cn; 


— SI -IS! CM 


CO— 


glg 




— > 


I *-H*— 'w*— • N ► 


uc 1 — ► 


— ►'^ ► 


►r>J 


— * — \ 


Q.JS 


X | 




— 


I <x<2: 


aC^- 


— Q — 


O0<w 


1 • — 


1 am- 


ca-w\— 


X «f — 1 — 




y! -0- — i*s.— • 


aJX 


x — >o — * x uj i — — — ►< x 


— CL — < 


— o<— < <£ 


<r\j2: 


— — ^ — -f 



jl.oc 
e ^ 

— h- 3 
I— *“ 

O—^O I 
X-IJJ021 
CCCLX^O II 

ox— ■— 2: 

'./)-< 0*sj;r 



2T \ + X I • — r\j\Xtt j ^\q + | r\j\ x 

•• + — — • — — •* • — . w— • C 9 | — * «- 

— i*C^ — r\j*: | oavNj<(\jpnrvj^Q-^:— *— i \ zi 

H — ~ >- || w || || || || h w i| || ^ n ~(\j || || || 

^xo— — .cy; — :c— x> i — — co co c\j — 
II l| II |l (\JcO<*0<\,rO |l 'NJ II CO II X 
r\j — - — * — - m — ► 



l -2 
— — 1 

— < II vC II II II II 
.O | — ~ — 

fM ii ii I — ii ii 

X 'O 1 — — 0- • ► ► * ► 



«r\jr\jr\jc\j ^y^yCCSXXXX II X:*C— i _i ^ r\j rO 4“ O 

—— .^^^wo-ww-ww n ii ——a ii ~ -x 

qcclxjq-c 4 <* < <j < ‘Oo<rti<32T r <<r'NZQ^NNJcauououj 



cm 



103 



L 



oooooooooooooooooo 



g 

8 

g 

0 

CfcJ 

£ 

H 

g 

1 

B 

W 

§2 

8 

a 



co 



I 

> 

I 



'i 



* 

i 




oooooooooooooooooo 

• •#••••#••#••••••• 

oooooooooooooooooo 



O' O' 0 s O' O'* O' O' ^ O' O' ^ O' O' ^ ^ ^ ^ ^ 

oooooooooooooooooo 

• •#•••••••#•••••#• 

oooooooooooooooooo 

O r^corornrnrncncncorocncorocncnfnfncn 

o 

pH OOOCOOCO^»tfl0^nOAvOONO 

<VI • •••••••••#••••••• 

in oo^tr-oooc7'(\jc\jsom(\j(\istm(\jcoc\j 
>t in>tO'cca'r^rof\i<\ja'a'Ovjc\irno'Ou^so 

iH 

o co in o fH O' cm cm ao oo o O' in cm cm o 
fnmcco s O'f s *>o»MHOo^oomninoin 

o • ••••• 

«j>t cncnrocMCMCNCMcnrococncncncncMcnr"^ 

~0' 

xo 

k OKCMmoinino'O'CMmin^coooom 

mcM *CvO*tcn^cMpHO'0'pH*HoococMKCD>fr 

• oooCst>t>toooocnfnrsjrsi^Ntf\j(M^oo 

<S)st • •«##••••••••••#•• 

i/) H CC^hhh^OOOOOOOOOOOOO 

UJ — 1 

JOO 

Z*’0fviin-<ino0'0--<o<\j<\ir-'0in>ooo#-ir- 

*-»ccoKKcn'OfMinin>$-*frCMCM>fr>a'KCMO'pHK 

< ^ o O' O' oa oo co in in m fn o^i ni cn ro nj (\i #-i o 

Kino # # 

l/>0 #000000000000000000 

• H 

oo 

m of\jooom-^ogor-i»-iap-i^K-^oo-^ 
** ^-in^co>t^^'0>oininKKKOooorsj 
000'0'rororo>f KKin^rorsj^ 

ooo 

«0'OMMN<NJ(VJH^OOOOOOOOOOO 

J'O^ 

Qp>H h 

oo ao •ouoouoo^Hoooooooino 

• o K oo >o in ro cm f-H cmt' <n m >r >y O' k o O' o 
ZO mcOKO'K>r >r '0'OrsirvJCTO^-^OOC7' 



cl KKr\jc\jro^)vonjrvjp-H^r-K^aoo<NJo 

H H p^p^p^ 

ocmo 

<<7'ooooooooooooinmooooo 
JOOOHmpnfn^'tmininin't^^'OOOO 
***0 o CO O' O' O' O' O' O' O' O' O' O' O' O' O' O' O' o o 

oo • ••••••••••♦••••••• 

2 *IMOOOOOOOOOOOQOOOOHp< 
<0-4 

*“* OOOOOOOOOOOOOOOOOO 
QC OOOOOOOOOOOOOOOOOO 
KOsJ • 

O' <\jc\jr\jinininininininininmmininoo 

#*CM H r « 

>rnj 

HOOOOOOOOOOOO<M(\IOJCVJMNN 
• O •OOOOOOOOOOO'OO'O^vOC'O 

O »OOOOOOOOOOOQ^hhhh«p4 

20 oooooooooooinininininnm 

>r ^>r^r<NrgcMCM04cv<\j 

at 

<oo 

«jr-o*MCMCMCMCMCMCM<MCMCMCMininininmifMn 
O^orifncnrnrnfnrnfnfnfnfriKKKKKKK. 
OOM'OfMCMCMCMCMCMCMCMfMCMCMKKKKKKK 
4*0 • • 



• —4 

o ^cMpris^m>CKcoc'o^cMco^rin'<jKop 

• ! < * ; i ,h(pIHHHp#pKiH # 



104 



INITIAL DISTRIBUTION LIST 



No. Copies 



1. Defense Documentation Center 20 

Cameron Station 

Alexandria, Virginia 22314 

2. Library 2 

Naval Postgraduate School 

Monterey, Calif. 93940 

3. Naval Ship Systems Command 1 

Department of the Navy 

Washington, D. C. 20360 

4. Professor P. F. Pucci 10 



Mechanical Engineering Department 
Naval Postgraduate School 
Monterey, Calif. 93940 

5. LCDR Henry J. Trost, Jr. 1 

Landing Ship Flotilla ONE 

FPO 

San Francisco, Calif. 96650 

6. Chairman 1 

Mechanical Engineering Department 

Naval Postgraduate School 
Monterey, Calif. 93940 



105 



Unclassified 

Security Classification 



DOCUMENT CONTROL DATA - R&D 

(Security ciaeellleatlon of title, body oi abetract and Indexing annotation muet be entered when the overall report la claeeUled) 



1. ORIGINATIN G ACTIVITY (Corporate author) 

Nsval Postgraduate School 
Monterey, California 93940 



la. REPORT SECURITY CLASSIFICATION 



26. GROUP 



3. REPORT TITLE 

The Use of the Analog Computer with the Single Blow Transient Testing Technique 
for Compact Heat Exchanger Surfaces 



4. DESCRIPTIVE NOTES (Type ot report and tnclvelve datee) 



5. AUTHORfSJ (Laet name, (l ret name, initial) 

TROST, Henry John Jr. 



f • REPORT DATE 



September 1967 



7a. TOTAL NO. OF 

104 



76. NO. OF REFS 

17 



• «. CONTRACT OR ORANT NO. 



9* ORIGINATOR** REPORT NUMBER^ 



6. PROJECT NO. 



N/A 



d. 



N/A 



§6. gTHER REPORT NOfSJ (Any other numbere that may be aaelfrod 



10. A VA IL ABILITY/LIMITATION NOTICES 




iy wLth approval of ISiaval Post* 



11. SUPPLEMENTARY NOTES 



12. SPONSORING MILITARY ACTIVITY 



13 ABSTRACT 

The purpose of this investigation was to determine the feasibility of 
using an analog computer to obtain the time derivative of the temperature 
response of a compact heat exchanger surface subjected to a step change in 
incoming fluid temperature; and to investigate the effect of the ratio of 
flow length to hydraulic diameter (L/D^) on the heat transfer and flow 
friction characteristics of compact heat exchanger surfaces. 

The method of maximum slope developed by Locke and modified by Howard 
to include conduction parameter was used to determine the heat transfer in- 
formation included herein. 

The results show that an analog computer can be a useful tool to aid 
in the collection and reduction of data. Results in the L/D investiga- 
tion were generally inconclusive and bear further investigation. 



DD 



FORM 

I JAN 64 



1473 



107 



Unclassif led 



Security Classification 



Unclassified 

- Security Classification 



key wo RDS 



ROLE 



ROLE 



Compact Heat Exchangers 
Flow Length to Hydraulic Diameter Ratio 
Maximum Slope by Analog Computer 
Differentiation on Analog Computer 



4 

i 



DD , F ° R :..1473 <back) 



S/N 0101 -807-6821 



Unclassified 



108 



Security Classification 



A- 3 1 409 



